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ABSTRACT 


An  investigation  of  steady  and  unsteady  flow  phenomena  in  centrifugal  compressors 
has  been  performed.  The  effect  of  vaned  diffuser  geometry  on  the  compressor  unsteady 
aerodynamics  has  been  considered,  with  particular  emphasis  on  the  diffuser  vane  unsteady 
loading  generated  by  the  impeller  circumferentially  nonuniform  flow  field.  A  series  of 
experiments  was  performed  in  the  Purdue  Research  Centrifugal  Compressor  Facility  to 
quantify  the  compres.»r  performance,  impeller  blade  and  diffuser  vane  steady  surface  static 
pressure,  vaneless  diffuser  steady  and  unsteady  velocity  field,  diffuser  vane  unsteady 
surface  static  pressure,  as  well  as  surge  and  rotating  stall  occurrence.  These  measurements 
were  made  for  various  flow  rates,  number  of  diffuser  vanes,  diffuser  vane  leading  edge- 
impeller  exit  radius  ratios,  diffuser  vane  stagger  angles,  and  nonuniform  circumferential 
vane  spacing. 

In  conjunction  with  the  above  experiments,  theoretical  predictions  of  the  unsteady 
vaneless  diffuser  wake  velocities  and  the  unsteady  diffuser  vane  loading  were  developed. 
The  linearized  Euler  equations  were  solved  to  predict  the  wake  behavior,  with  a  conformal 
transformation  applied  to  existing  axial  flow  cascade  theory  for  unsteady  diffuser  vane 
loading  predictions. 

It  was  found  that  the  compressor  performance  could  be  improved  through  use  of  a 
vaned  diffuser,  with  the  greatest  performance  improvement  with  zero  steady  loading  on  the 
vanes.  The  unsteady  wake  velocity  in  the  radial  diffuser  was  seen  to  decrease  rapidly  with 
increasing  radius.  The  decrease  was  greater  than  predicted  by  the  theory,  with  the  unsteady 
circumferential  pressure  gradient  and  the  phase  angle  between  the  radial  and  tangential 
unsteady  velocities  determined  to  be  important  parameters. 

The  diffuser  vane  unsteady  loading  was  determined  to  be  a  strong  function  of  flow 
rate  and  number  of  diffuser  vanes,  with  the  diffuser  vane  leading  edge-impeller  exit  radius 
ratio  important  for  smaller  numbers  of  vanes.  The  correlation  with  theory  was  poor  with  a 
high  diffuser  vane  row  solidity,  although  fair  correlation  was  seen  with  moderate  solidity. 

Finally,  the  vaned  diffuser  row  had  a  significant  impact  on  surge  and  rotating  stall 
margins  and  behavior.  It  was  found  that  notable  improvement  could  be  made  by  proper 
adjustment  of  the  diffuser  vane  row  geometry. 
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CHAPTER  1 
INTRODUCTION 


Centrifugal  compressors  have  found  wide  use  in  industrial,  aerospace  and  ground 
based  transportation  applications.  Important  systems  include  small  gas  turbine  engine 
compressors,  piston  engine  turbochrrgers,  compressed  air  and  refrigeration  systems,  and 
turbopumps  for  liquid  fueled  rocket  engines.  The  earliest  jet  engines  employed  centrifugal 
compressors.  However,  axial  compressors  became  favored  due  to  their  smaller  cross 
sectional  area  which  facilitated  use  in  high-thrust  jet  engines.  Centrifugal  compressors  do 
have  several  distinct  advantages,  including  their  relative  ease  of  manufacture  and  much 
higher  pressure  rise  per  stage  as  compared  with  axial  machines.  These  advantages  have 
begun  to  be  exploited  in  the  past  two  decades,  prompting  the  need  for  research  to  improve 
centrifugal  compressor  performance. 

A  centrifugal  compressor  is  a  rotating  machine  which  increases  the  stagnation  enthalpy 
of  a  flowing  fluid,  with  the  throughflow  mainly  in  a  radial  direction.  It  is  this  radial 
throughflow  which  allows  a  centrifugal  compressor  to  attain  a  very  high  pressure  rise  per 
stage.  The  impeller  and  radial  diffuser  comprise  a  centrifugal  compressor  stage,  with  inlet 
guide  vanes  sometimes  utilized  to  impart  prewhirl  to  the  incoming  flow.  The  radial  diffuser 
is  generally  vaneless.  However,  many  modem  compressor  stages  employ  vaned  diffusers 
to  increase  efficiency. 

Two  typical  centrifugal  compressor  stages  are  depicted  in  Figure  1.1.  The  fluid  enters 
the  impeller  axially  at  the  eye  or  inlet  Ideally,  the  forward  portion  of  the  impeller,  termed 
the  inducer,  allows  the  flow  to  enter  the  impeller  smoothly  and  directs  the  flow  radially. 
The  fluid  static  pressure  is  continuously  increased  throughout  the  radial  portion  of  the 
impeller  due  to  centripetal  acceleration.  A  static  pressure  gradient  exists  across  the  impeller 
channel  resulting  from  the  coriolis  acceleration  of  the  fluid.  This  gradient  produces  a 
pressure  difference  across  each  impeller  blade  which  exerts  torque  on  the  impeller.  The 
flow  leaving  the  impeller  attains  a  greater  velocity  than  the  inlet  flow  due  to  the  added  wheel 
speed.  The  radial  diffuser  is  designed  to  receive  this  flow  and  to  attain  the  largest  possible 
pressure  recovery. 
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The  actual  flow  through  a  centrifugal  compressor  stage  is  highly  complex,  as  depicted 
schematically  in  Figure  1.2  (Wood,  et  al.,  1983).  The  flow  through  the  impeller  passages 
is  complicated  by  several  interacting  phenomena.  The  static  pressure  and  velocity  are 
circumferentially  nonuniform  at  the  impeller  exit  due  to  the  potential  effects  of  flow  through 
the  impeller  passage  and  viscous  effects  manifested  in  trailing  edge  wakes  and  flow 
separation.  This  circumferentially  nonuniform  flow  is  unsteady  to  the  stationary  diffuser, 
which  is  expected  to  attain  the  largest  possible  pressure  recovery.  The  potential  disturbance 
of  the  diffuser  vanes  is  also  felt  upstream,  imposing  an  unsteady  backpressure  on  the 
rotating  impeller.  In  addition  to  these  periodic  disturbances  with  the  rotational  or  blade 
passing  frequency,  instabilities  can  develop  in  centrifugal  compression  systems  in  the  form 
of  rotating  stall  and  surge.  All  of  these  unsteady  interactions  and  disturbances  can 
significantly  impact  the  compre::or  performance  and  produce  aerodynamically  induced 
stress  and  vibration  which  may  result  in  structural  failure. 

It  has  therefore  become  necessary  to  attain  a  greater  understanding  of  the  unsteady  flow 
phenomena  in  centrifugal  compressors.  Since  many  modem  centrifugal  compressors 
employ  vaned  diffusers,  experimental  investigations  in  conjunction  with  analyses  are 
needed  to  produce  the  most  efficient  designs  and  minimize  undesirable  unsteady  flow 
phenomena. 


U  Operation  p£.Csntrifugal  Impellers 

The  flow  enters  the  forward  portion  of  the  impeller,  which  is  termed  the  inducer.  It  is 
designed  to  allow  the  flow  to  enter  the  impeller  smoothly  and  to  turn  the  flow  to  a  radial 
direction  with  minimal  loss.  The  inducer  inlet  blades  are  typically  given  some  angle  with 
respect  to  axial  to  match  the  relative  angle  of  the  inlet  flow  over  a  certain  range  of  operation. 
Most  of  the  energy  transfer  to  the  air  occurs  in  the  radial  portion  of  the  impeller.  The 
blades  at  the  impeller  exit  typically  are  radially  oriented  or  given  some  angle  from  radial 
opposite  the  impeller  rotation.  The  motivation  for  radially  oriented  blades  includes  ease  of 
manufacture  and  reduced  stress  in  high  speed  machines.  Backswept  blades  are  used  to 
stabilize  the  impeller  flow  and  reduce  the  velocity  into  the  diffuser  which  aids  in  pressure 
recovery.  This  lowers  the  power  output  of  the  stage  but  can  increase  the  efficiency. 
Forward  swept  blades  are  rarely  used. 

In  a  centrifugal  impeller,  as  previously  mentioned,  the  flow  in  the  impeller  channel  will 
be  nonuniform,  producing  a  temporal  variation  in  both  the  pressure  and  velocity  at  the 
impeller  exit.  This  causes  the  flow  in  the  radial  diffuser  to  be  unsteady.  In  fact,  Dean 
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(1959)  shows  that  it  is  impossible  for  a  turbomachine  to  work  on  a  flowing  fluid  by  other 
than  viscous  forces  without  time  variation  in  the  flow.  It  is  therefore  impossible  for  a 
turbomachine  to  operate  without  unsteady  flow  phenomena,  with  these  phenomena 
generally  increasing  in  magnitude  as  the  work  output  of  the  machine  increases. 


1.2  Basic  Theory  of  Radial  Diffusers 

The  radial  diffuser  is  a  critical  component  in  the  centrifugal  compressor  stage.  Its 
purpose  is  to  reduce  the  velocity  from  the  impeller  exit  and  to  increase  the  pressure  as 
efficiendy  as  possible.  Vaned  diffusers  are  typically  employed  to  increase  efficiency. 

It  is  helpful  both  in  the  study  of  vaneless  and  vaned  diffusers  to  develop  a  basic 
analysis  for  flow  in  a  vaneless  diffuser.  The  continuity  equation  and  the  radial  and  angular 
momentum  equations  for  axisymmetric  incompressible  frictionless  flow  in  a  radial  parallel 
wall  diffuser  are 


fiWUr)  =  0 

rdr  1.10 

1.11 

1.12 

where  the  flow  is  assumed  to  be  steady  for  the  simplified  analysis  and  the  tangential 
pressure  gradient  is  assumed  to  be  zero. 

Integrating  Equations  1.10  and  1.12  yields 


r  dr  r  Po  dr 

Uf§U£  +  UrUsUo 

OT  r 


rUr  =  riUrl 

1.13 

rUe  =  riUei 

1.14 

where  the  subscript  1  refers  to  the  diffuser  inlet 

Combining  Equations  1.11, 1.13  and  1.14  and  integrating  yields 
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IflPoUf  1.15 

where  uf  =  U§|  +  U^.  Thus,  the  pressure  decreases  with  the  inverse  of  the  radius 
squared. 

In  an  incompressible  flow,  both  the  radial  and  the  tangential  velocity  components 
decrease  as  1/r.  Thus,  from  continuity  and  conservation  of  angular  momentum,  both 
velocity  components  are  reduced  in  the  vaneless  diffuser  due  to  the  increase  in  radius. 

It  is  useful  in  the  understanding  of  diffuser  flows  to  obtain  the  equation  for  a 
streamline.  The  equation  of  a  streamline  for  an  incompressible  flow  is  given  by 

lilL=UL 

rd0  u0  1.16 

Assuming  steady,  frictionless,  axisymmetric  flow.  Equations  1.13  and  1.14  can  be  used 
along  with  Equation  1.16  yielding 


ljdL=  U-iL 

r  d0  uei  1.17 

where  the  subscript  1  again  refers  to  the  diffuser  inlet 
Integrating  Equation  1.17  from  diffuser  inlet  yields 

ln-L  =  (e-0i)/tanaL  ug 

where  the  angle,  measured  from  the  radial  direction,  oil  =  tan'^Uei/Uri)  and  the  subscript 
1  refers  to  the  diffuser  inlet 

Equation  1.18  is  the  equation  of  a  logarithmic  spiral,  which  forms  an  angle  at  with  any 
radial  line,  where  gcl  is  termed  the  logarithmic  spiral  angle.  Thus,  the  incompressible  flow 
in  a  parallel  wall  radial  diffuser  follows  logarithmic  spiral  streamlines  whereas  the  flow  in 
an  axial  machine  tends  to  follow  a  line  with  a  constant  angle  from  axial.  This  fact  must  be 
considered  in  the  analysis  and  design  of  radial  turbomachine  components. 
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1.3  Literature  Review 

Much  attention  in  the  literature  has  been  given  to  centrifugal  impeller  and  diffuser 
flowfields.  The  impeller  exit  flow,  depicted  in  Figure  1.3,  has  typically  been  described  as 
a  jet- wake  type  of  flowfield.  Because  of  boundary  layer  growth  and  flow  separation  in  the 
impeller,  the  exit  velocity  profile  is  circumferentially  nonuniform.  This  type  of 
circumferentially  nonuniform  flow  has  been  discussed  by  several  authors.  Dean  and  Senoo 
(1960)  formulated  a  mathematical  model  of  this  jet- wake  impeller  flow.  They  concluded 
that  a  reversible  work  interaction  in  circumferentially  nonuniform  flows  between  the  jet  and 
wake  causes  the  velocity  deficit  to  decrease  rapidly  within  a  small  radial  distance  from  the 
impeller  exit  They  also  discovered  that  the  unsteady  tangential  velocity  component  in  the 
radial  diffuser  can  increase  before  its  eventual  decay.  Eckardt  (1975)  performed 
measurements  with  high  frequency  instrumentation  in  the  jet-wake  flow  of  a  low  pressure 
ratio  centrifugal  impeller.  He  found  a  significant  jet-wake  pattern  existed  and  that  the  radial 
extension  of  the  nonuniform  circumferential  velocity  profile  decayed  more  slowly  than 
predicted  by  mathematical  models  such  as  the  one  formulated  by  Dean  and  Senoo. 

Krain  (1987)  measured  the  velocity  field  within  a  newly  designed  impeller  using  an 
L2F  system.  Velocity  profile  nonuniformities  were  found  within  the  impeller  passages  but 
a  relatively  uniform  velocity  profile  at  the  impeller  exit  which  differed  significantly  from  the 
typical  jet  wake  pattern.  Fagan  and  Fleeter  (1989)  measured  the  flowfield  in  the  passages 
of  a  mixed-flow,  low  pressure  ratio  centrifugal  impeller.  The  traditional  jet-wake  exit 
flowfield  pattern  was  observed  at  the  design  point  and  near  stall.  From  these 
investigations,  it  appears  that  although  the  typical  impeller  exit  flow  is  nonuniform 
circumferentially,  the  actual  flow  pattern  is  dependent  on  the  impeller  design,  with  even 
near  uniform  profiles  possible  if  separation  in  the  impeller  is  controlled.  The  many 
interacting  phenomena  in  impeller  flows  are  not  well  understood,  however,  and  it  does  not 
seem  reasonable  to  expect  either  complete  elimination  of  separation  or  circumferentially 
nonuniform  flows  in  radial  impellers. 

Vaned  diffusers  serve  as  flow  guides  and  also  incorporate  turning  of  the  flow  so  as  to 
reduce  the  tangential  velocity  component  at  a  greater  rate  than  in  a  vaneless  diffuser,  thus 
reducing  the  size  of  the  diffuser.  Vaned  diffusers  have  typically  been  of  two 
configurations:  cambered  vane  diffusers  and  wedge  or  channel  diffusers.  Representative 
examples  of  these  diffusers  are  shown  in  Figure  1.1.  The  cambered  vane  diffuser  has 
developed  from  analogy  to  airfoil  or  turning  vane  cascades.  Cambered  vane  diffusers  with 
logarithmic  spiral  profiles  and  little  or  no  flow  turning  have  been  found  to  increase  stage 
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efficiency  (Smith,  1970).  This  is  attributed  to  their  ability  to  inhibit  secondary  flows  and 
separation.  Thus,  use  of  vaned  diffusers  simply  as  flow  guides  can  increase  performance. 

The  second  major  type  of  vaned  diffuser  is  the  wedge  or  channel  diffuser.  The  purpose 
of  the  channel  diffuser  is  to  remove  the  tangential  velocity  component  and  to  limit  the  area 
increase  of  the  diffuser.  Channel  diffusers  are  generally  designed  using  some  divergence 
angle  obtained  from  two-dimensional  straight  diffuser  theory.  Channel  diffusers  have  been 
given  much  attention  in  the  literature  and  have  shown  good  performance. 

Clements  and  Artt  (1986)  and  Rodgers  (1982)  investigated  the  performance  of  channel 
type  diffusers  considering  the  diffusion  capabilities  of  different  geometries.  Clements  and 
Artt  conclude  that  the  performance  of  channel  diffusers  is  superior  to  that  of  cambered  vane 
diffusers.  Bammert,  Jansen  and  Rautenberg  (1983)  considered  both  types  of  diffusers 
along  with  a  twisted  vane  diffuser.  Their  data  does  not  support  the  conclusion  of  Clements 
and  Artt,  with  cambered  vane  diffuser  stages  having  somewhat  higher  efficiency. 
Yoshinaga  et  al.  (1980)  investigated  the  performance  of  sixteen  cambered  vane  diffusers 
with  moderate  to  high  diffusion  ratio.  They  compared  their  experimental  results  with 
potential  flow  solutions  using  conformal  mapping  and  found  fair  to  good  correlation  with 
the  theory. 

The  effect  of  the  radius  ratio  of  the  vaneless  space  in  vaned  radial  diffusers  has  been 
discussed  by  several  authors.  Optimum  values  given  range  from  1.1  to  1.15  (Rodgers, 
1982;  Bammert  et  al.,  1983;  Fisher  and  Inoue,  1981;  Inoue  and  Cumpsty,  1984).  It  is 
generally  concluded  in  the  literature  that  vaned  diffusers  become  less  effective  at  larger 
radius  ratios.  Therefore,  it  is  desired  to  position  the  diffuser  vanes  as  close  as  possible  to 
the  impeller  while  avoiding  large  unsteady  interactions  and  reduced  performance  caused  by 
these  interactions. 

Very  few  investigations  have  considered  the  unsteady  flow  phenomena  such  as  the 
interactions  between  the  impeller  and  the  vaned  diffuser  or  unsteady  diffuser  vane 
performance.  Fisher  and  Inoue  (1981)  made  instantaneous  measurements  in  the  vaneless 
space  of  a  low  speed  compressor  stage.  However,  only  time  averaged  data  were  compared 
to  steady  flow  predictions  from  conformal  mapping  and  Martensen  singularity  theory 
(Gostelow,  1984),  with  good  correlation  obtained.  Although  several  authors  have 
concluded  that  the  vaned  diffuser  had  little  effect  on  the  overall  impeller  performance, 
Fisher  and  Inoue  concluded  that  the  vaned  diffuser  significantly  affects  the  impeller 
d'  charge  and  vaneless  space  flowfield.  Inoue  and  Cumpsty  (1984)  investigated  impeller 
discharge  flow  in  vaneless  and  vaned  diffusers.  The  vaned  diffusers  employed  10, 20  and 
30  blades  and  were  operated  at  diffuser  entrance/impeller  tip  radius  ratios  of  1.04, 1.1  and 
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1.2.  They  found  circumferentially  nonuniform  velocity  profiles  in  the  vaneless  diffuser  at 
large  radii,  but  these  nonuniformities  did  not  propagate  very  far  into  the  vaned  diffuser. 
Arndt  et  al.  (1988)  investigated  the  interaction  between  impeller  blade  and  diffuser  vanes  in 
a  diffuser  pump.  Their  measurements  were  made  at  various  flow  rates,  impeller  speeds, 
and  diffuser  vane  leading  edge/impeller  exit  radius  rados.  They  found  increasing  the  radius 
ratio  strongly  decreased  the  unsteady  vane  loading. 

In  the  above  unsteady  diffuser  flow  investigations,  only  Arndt  measured  the  unsteady 
vane  loading,  with  no  comparison  with  theory.  In  fact,  no  general  theory  for  unsteady 
vane  loading  due  to  impeller-diffuser  interactions  is  known.  Therefore,  it  is  important  to 
extend  the  experimental  data  base  of  unsteady  impeller-diffuser  interactions,  with  detailed 
flowfield  measurements,  as  well  as  developing  a  theory  to  predict  these  interactions. 

Greitzer  (1981)  discussed  the  types  of  instabilities  that  arise  in  pumping  and 
compression  systems.  He  showed  analytically  that  system  dynamic  instabilities  in  the  form 
of  surge  can  develop  in  compression  systems  when  the  slope  of  the  pressure  characteristic 
becomes  positive.  He  also  considered  rotating  stall  in  centrifugal  compressors  and  showed 
that  rotating  stall  inception  can  occur  before  surge,  although  the  mechanisms  for  rotating 
stall  inception  are  not  well  understood.  Van  den  Braembussche  (1984)  considered  rotating 
stall  and  surge  in  centrifugal  compressors.  He  believes  that  rotating  stall  occurs  as  a  self- 
sustaining  oscillation  in  which  the  inviscid  flow  field  reacts  to  a  viscous  perturbation, 
typically  stalling  of  the  impeller  or  vaned  diffuser. 

Very  little  of  the  literature  concerning  rotating  stall  and  surge  in  centrifugal  compressors 
considers  diffuser  geometry.  It  is  believed  that  diffuser  geometry  can  have  an  important 
effect  on  this  behavior.  Therefore,  it  is  desirable  to  investigate  the  effect  of  diffuser 
geometry  on  surge  and  rotating  stall. 

Although  the  literature  concerning  centrifugal  compressor  diffusers  is  quite  extensive, 
the  emphasis  has  been  on  steady  flow  performance,  with  very  few  comparisons  between 
predictions  and  data.  Examples  of  these  comparisons  include  Yoshinaga  (1980),  and 
Fisher  and  Inoue  (1981).  Faulders  (1956)  discusses  the  application  of  potential  flow 
theory  and  conformal  mapping  to  existing  steady  cascade  data  to  predict  diffuser 
performance.  Smith  (1970)  calls  fe*  a  diffuser  technology  based  on  conformal 
transformation  without  regard  to  channel  or  cambered  vane  types.  That  a  coherent 
technology  and  method  of  presentation  of  data  is  lacking  is  evident  from  the  various  and 
sometimes  contradictory  conclusions  drawn  by  different  authors.  Although  Faulders  paper 
was  published  in  1956,  very  few  authors  have  attempted  to  correlate  data  using  his  work. 
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Much  of  the  variation  in  the  published  data  is  probably  due  to  the  dependence  of  the  data  on 
the  particular  research  facility. 

Bryan  and  Fleeter  (1987)  experimentally  investigated  the  effects  of  inlet  prewhirl  on  the 
impeller  flowfield  and  overall  performance  of  the  Purdue  Research  Centrifugal 
Compressor.  This  is  a  large  scale,  low  speed  machine  which  was  originally  designed  and 
operated  by  Vavra  (1955).  The  facility  was  established  at  Purdue  by  Bryan  and  Fleeter  and 
is  used  in  the  current  investigation.  It  was  shown  that  inlet  prewhirl  significantly  affected 
the  impeller  flow,  especially  in  the  inducer  region  and  that  performance  could  be  optimized 
by  varying  inlet  prewhirl.  It  was  also  discovered  that  very  little  effect  of  the  prewhirl 
propagated  into  the  impeller  discharge/diffuser  region  where  the  flowfield  was  almost 
solely  determined  by  flow  rate. 

Thus  in  summary,  considerable  effort  has  been  made  to  determine  the  performance  of 
centrifugal  compressor  stages  with  emphasis  on  the  radial  diffuser.  Investigations  of 
unsteady  phenomena  and  the  effect  of  diffuser  geometry  on  such  phenomena,  however, 
have  bten  given  relatively  little  attention,  with  few  or  no  theoretical  predictions.  Therefore, 
it  is  imperative  that  these  unsteady  phenomena  be  investigated  and  the  important  machine 
parameters  affecting  them  be  determined. 

1.4  Research  Objective  and  Technical  Approach 

The  overall  objective  of  this  investigation  is  to  quantify  the  dependence  of  significant 
steady  and  unsteady  flow  phenomena  inherent  in  centrifugal  compressors  on  vaned  diffuser 
geometry,  with  emphasis  on  the  radial  vaned  diffuser  unsteady  aerodynamics  generated  by 
the  impeller-diffuser  interaction.  Specific  objectives  include: 

•  Determination  of  Centrifugal  Compressor  Performance 

•  Impeller  and  Diffuser  Vane  Steady  Loading  Measurement 

•  Diffuser  Unsteady  Flowfield  Measurement 

•  Diffuser  Vane  Unsteady  Loading  Measurement 

•  Determination  of  Surge  and  Rotating  Stall  Occurrence  and  Characteristics 

•  Development  of  Theory  to  Determine  Radial  Vaneless  Diffuser 
Wake  Behavior  tnd  Unsteady  Diffuser  Vane  Loading 

The  above  experimental  objectives  were  accomplished  by  performing  a  series  of 
experiments  in  the  Purdue  Centrifugal  Compressor  Facility.  Steady  flowfield 
measurements  were  made  using  pressure  probes  throughout  the  machine  to  determine 
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performance.  Steady  surface  static  pressure  measurements  were  made  on  the  rotating 
impeller  and  vaned  diffuser  to  determine  loading.  The  time-variant  flowfield  in  the  radial 
diffuser  were  measured  using  an  LDV  system.  High  response  pressure  transducers  were 
utilized  to  measure  the  unsteady  diffuser  vane  surface  static  pressure  as  well  as  static 
pressure  variations  at  the  impeller  exit  for  determination  of  surge  and  rotating  stall. 

The  above  experiments  were  performed  for  various  diffuser  geometries  and  flow  rates. 
In  particular,  the  number  of  diffuser  vanes,  the  diffuser  vane  stagger,  the  diffuser  vane 
leading  edge/impeller  exit  radius  ratio  and  the  circumferential  spacing  were  varied 
independently  along  with  the  mass  flow  rate  through  the  machine.  The  compressor 
performance  and  detailed  flowfield  were  determined  at  each  machine  operating  condition. 

Analyses  were  developed  using  linearized  theory  to  predict  the  unsteady  wake  behavior 
in  the  radial  vaneless  diffuser.  This  was  incorporated  into  existing  unsteady  cascade  theory 
using  conformal  mapping  to  predict  the  unsteady  loading  on  the  vaned  diffuser  due  to 
wakes  shed  from  the  rotating  impeller. 
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Figure  1.1  Typical  Centrifugal  Compressor  Stages 
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CHAPTER 2 

THE  PURDUE  RESEARCH  CENTRIFUGAL  COMPRESSOR 

2.1  General 

The  Purdue  Research  Centrifugal  Compressor  experimentally  models  the  fundamental 
aerodynamic  phenomena  inherent  in  centrifugal  compressor  stages  including  the 
circumferentially  nonuniform  impeller  exit  flow,  diffuser  vane  incidence  and  radius  ratio, 
and  unsteady  impeller  diffuser  interactions.  The  compressor  is  belt  driven  using  a  jackshaft 
and  pulley  arrangement  by  a  29.8  kW  (40  hp)  induction  motor  and  is  operated  at  1790  rpm. 
The  facility  is  shown  in  Figure  2. 1  and  depicted  schematically  in  Figure  2.2. 

The  general  operation  of  the  machine  can  be  described  by  considering  the  flow  path 
shown  in  Figure  2.3.  Air  is  drawn  into  the  machine  through  the  bellmouth  inlet  section 
which  converges  to  an  axial  direction  at  the  impeller  eye.  The  air  enters  the  mixed  flow 
impeller  which  changes  the  axial  inlet  flow  to  a  nearly  radial  flow  at  the  outlet  and  also 
imparts  a  tangential  component  of  velocity  to  the  air.  The  center  of  the  rotor  is  hollow, 
permitting  the  installation  of  instrumentation  to  transmit  pressure  data  from  the  surfaces  of 
the  rotating  blades  to  the  data  acquisition  and  analysis  system  located  outside  of  the 
macnine  in  the  stationary  frame  of  reference. 

On  leaving  the  impeller,  the  air  enters  the  radial  diffuser  whose  walls  curve  smoothly  to 
a  parallel  configuration.  The  air  then  passes  through  a  row  of  stationary  diffuser  vanes 
which  are  located  at  the  upstream  portion  of  the  parallel  wall  diffuser  and  is  discharged  into 
a  large  annular  plenum  chamber.  From  this  plenum,  the  air  passes  into  a  straight  discharge 
pipe  which  contains  a  sharp  edged  orifice  for  flow  rate  measurement  and  a  throttle  valve  to 
regulate  the  flow  rate. 


2?,  Impeller, and  firivaAssfimbly 

The  impeller  and  drive  assembly  are  shown  in  Figure  2.4.  The  impeller  consists  of  23 
aluminum  blades  cast  integrally  with  the  hub.  The  tip  is  covered  with  a  plexiglass  shroud 
for  aerodynamic  purposes  and  optical  access.  A  mixed-flow  impeller  configuration  is  used 
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with  a  deloaded  blade  profile.  In  a  mixed-flow  impeller,  the  flow  channels  never  turn  to  a 
completely  radial  direction.  Therefore,  the  flow  maintains  an  appreciable  component  of 
velocity  in  the  axial  direction  during  its  entire  passage  through  the  impeller.  With  an 
impeller  of  this  type,  the  blade  need  not  lie  in  a  purely  meridional  plane,  but  may  follow  a 
spiral  path  along  the  hub  from  inlet  to  exit  The  use  of  a  spiral  type  of  flow  path  introduces 
a  great  deal  of  flexibility  into  the  aerodynamic  design  and  permits  use  of  a  deloaded  blade 
design.  In  a  general  way,  this  is  an  impeller  blade  on  which  the  pressure  side  is  concave 
over  the  first  part  of  its  length  and  essentially  convex  over  the  last  part.  This  reversal  of 
curvature  tends  to  reduce  the  aerodynamic  loading  near  the  blade  exit  and  increase  the 
loading  over  the  central  portion  of  the  blade. 

The  steady  pressure  distributions  in  the  impeller  blade  passages  are  measured  with 
static  pressure  taps  on  the  blade  surfaces.  There  are  32  such  taps,  with  two  redundant 
pairs,  located  along  three  streamlines  at  five  normal  locations.  These  are  shown  in  Figure 
2.5.  The  steady  pressures  on  the  impeller  are  measured  using  a  rotating  pressure 
multiplexer. 

The  shaft  and  motor  are  track  mounted  on  a  steel  base.  The  entire  impeller-drive 
assembly  can  be  retracted  to  perform  necessary  service.  The  impeller  assembly  utilizes 
four  cam  mounted  wheels  which  are  lowered  onto  the  track  upon  retracting  the  impeller. 
When  the  impeller  is  in  place,  the  wheels  are  raised  allowing  the  base  to  rest  in  a  v-groove 
in  the  track.  The  base  and  impeller  drive  assembly  are  then  bolted  to  the  track.  The 
impeller  is  fixed  in  the  stationary  passages  by  a  metal  seal  which  bolts  onto  the  rear  passage 
of  the  diffuser  section  downstream  of  the  impeller.  The  inner  part  of  this  seal  contains  the 
rear  labyrinth  seal. 


2.3  Inlet 

The  facility  inlet  allows  air  to  flow  smoothly  from  the  atmosphere  into  the  impeller. 
The  inlet  tip  wall  consists  of  a  plexiglass  bellmouth  of  identical  profile  to  the  impeller  blade 
tip.  The  inlet  hub  or  centerbody  consists  of  a  conical  section  smoothly  tapered  to  a  cylinder 
with  a  diameter  equal  to  that  of  the  impeller  hub.  The  inlet  has  been  developed  using  wall 
tufts  to  detect  separation  points  or  flow  nonuniformities.  Pressure  probe  data  show  the 
flow  to  be  ve:  y  uniform  at  the  impeller  face. 


15 


2.4  Diffuser 

The  radial  diffuser  consists  of  a  smoothly  curved  vaneless  space  followed  by  a  parallel 
wall  vaned  radial  diffuser.  The  diffuser  vanes  have  a  chord  length  of  16.5  cm  (6.5  in.), 
with  a  NACA  4312  airfoil  profile.  Figure  2.6.  The  wall  profile.  Figure  2.7,  facilitates  the 
close  spacing  of  the  diffuser  vane  leading  edge  and  impeller  trailing  edge.  This  profile  was 
obtained  by  fitting  a  cubic  spline  from  the  impeller  trailing  edge  to  the  parallel  wall  radial 
diffuser.  The  performance  of  the  diffuser  flow  path  was  evaluated  by  computing  the 
axisymmetric  potential  flow  and  applying  the  boundary  layer  code  CBTSL  with  a  Mangier 
transformation,  utilizing  impeller  exit  flowfield  data  obtained  by  Fagan  (1989)  as  a 
boundary  condition.  The  flow  in  'he  diffuser  was  smooth  and  nonseparated  over  the 
design  operating  range,  as  indicated  through  flow  visualization  using  wall  tufts.  The 
diffuser  tip  wall  contains  30  eccentric  blade  mounting  cams,  Figure  2.8,  which  allow  the 
diffuser  vane  leading  edge  radius  and  stagger  angle  to  be  varied  and  also  permit 
nonuniform  circumferential  vane  spacing. 

The  diffuser  walls  are  instrumented  with  a  field  of  static  pressure  taps  as  well  as  static 
taps  along  the  diffuser  vane  chord  at  mid-span.  A  three  head  cobra  probe  is  situated  at  the 
impeller  exit  to  measure  total  and  static  pressure  as  well  as  flow  angle.  Two  similar  probes 
are  located  downstream  of  the  diffuser  vane  row. 


2.5  Elsnunuod  Exhaust  Piping 

An  overhead  schematic  of  the  complete  facility  is  shown  in  Figure  2.9.  The  diffuser 
section  is  attached  to  the  plenum  into  which  the  flow  is  discharged.  The  plenum  is  toroidal 
with  a  square  cross  section  The  inner  diameter  is  1.69  m  (66.4  in.),  the  outer  diameter  is 
3.15  m  (123.9  in.),  and  the  width  is  74.2  cm  (29.2  in.).  The  plenum  has  a  circular  outlet 
with  a  70.0  cm  (24  in.)  i.d.  The  exhaust  piping  system  is  attached  to  the  plenum  exit.  As 
air  leaves  the  plenum,  it  enters  a  90°  elbow  with  an  inner  diameter  of  70.0  cm  (24.0  in.). 
The  elbow  is  followed  by  a  reducer  which  reduces  the  inner  diameter  to  40.6  cm  (16  in.). 
A  straight  section  of  40.6  cm  (16  in.)  inner  diameter  and  18  pipe  diameters  in  length 
follows  the  reducer.  A  flow  straightener  of  7.6  cm  (3  in.)  diameter  tubing  and  55.9  cm  (22 
in.)  length  is  located  2  pipe  diameters  downstream  of  the  reducer  to  remove  any  swirling  or 
radial  nonuniformities  in  the  flow.  A  sharp  edged  orifice  plate  of  30.48  cm  (12.0  in.) 
orifice  diameter  in  conjunction  with  1-diameter  and  1/2-diameter  pressure  taps  is  located  at 
the  end  of  the  straight  section  for  mass  flow  measurement.  A  straight  section  8  pipe 
diameters  in  length  is  located  behind  the  orifice  plate  and  contains  a  thermocouple  6  pipe 
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diameters  downstream  of  the  orifice.  A  butterfly  type  throttle  valve,  driven  by  a  gear 
motor,  is  located  at  the  exit  of  the  90°  elbow  to  control  the  compressor  flow  rate.  The 
motor  is  actuated  through  a  set  of  relays  under  control  of  the  data  acquisition  system. 


Figure  2. 1  Front  View  of  Compressor  Test  Facility 


Figure  2.2  Cutaway  of  Compressor  Test  Facility 
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Figure  2.3  Inlet,  Impeller,  and  Diffuser  Flow  Channel 
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Figure  2.7  Diffuser  Wail  Profile 
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Figure  2.9  Schematic  of  Centrifugal  Compressor  Facility 


CHAPTER  3 

DATA  ACQUISITION  AND  ANALYSIS 


In  the  current  investigation,  the  data  of  fundamental  interest  are:  1)  the  compressor 
steady  operating  conditions,  including  the  aerodynamic  performance  of  the  impeller  and 
diffuser;  2)  the  detailed  impeller  blade  and  diffuser  vane  steady  aerodynamic  loading 
distributions,  defined  by  the  surface  static  pressure  distributions;  3)  the  diffuser  mean  flow 
field;  4)  the  impeller  exit-diffuser  vane  inlet  mean  and  periodic  flow  field;  3)  the  unsteady 
pressure  distributions  on  the  diffuser  vane  surfaces;  and  6)  the  unsteady  static  pressure  at 
the  impeller  exit/diffuser  entrance. 

In  this  chapter,  the  instrumentation,  data  acquisition  system  and  analysis  techniques  to 
acquire  these  data  are  discussed.  Static  pressure  taps  on  the  impeller  blades  and  diffuser 
vanes  and  high  response  miniature  pressure  transducers  embedded  in  the  diffuser  vanes 
were  used  to  determine  the  detailed  steady  and  unsteady  aerodynamic  data.  The  diffuser 
steady  and  time-variant  flow  Held  information  was  acquired  with  direction  sensitive 
pressure  probes  and  laser  Doppler  velocimctery  (LDV). 


3.1  Steady  Data  Instrumentation 

The  centrifugal  compressor  steady  data  acquisition  system  is  described  in  detail  by 
Bryan  and  Fleeter  (1987).  The  time-averaged  pressures  on  the  stationary  components  are 
measured  using  6.9  kPa  (1  psi)  strain-gage  pressure  transducers  in  conjunction  with  a 
Scanivalve  pressure  multiplexing  system.  The  transducers  are  calibrated  using  a  water 
manometer  with  reference  pressure  and  vacuum  sources.  The  transducer  output  voltage  is 
measured  via  computer  using  a  Hewlett  Packard  HP3497  Data  Acquisition  Unit. 

The  impeller  is  instrumented  with  32  static  pressure  taps  on  the  rotating  blade  surfaces. 
The  static  pressures  are  measured  using  a  rotating  Scanivalve  pressure  multiplexer  and 
strain  gage  pressure  transducer.  The  rotating  Scanivalve  multiplexer  is  capable  of 
measuring  36  static  pressures  on  a  rotating  impeller.  The  Scanivalve  employs  a  pneumatic 
stepping  motor  and  optical  encoder  for  positioning.  The  pressure  from  the  rotating 
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multiplexer  assembly  is  exposed  to  a  stationary  pressure  transducer  by  way  of  a  rotating  o- 
ring  gland  seal.  The  Scanivalve  is  connected  to  the  impeller  through  a  manifold  and 
flexible  coupling.  The  entire  assembly  is'shown  in  Figure  3.1. 

The  rotating  Scanivalve  transducer  is  mounted  on  the  compressor  axis  while  the 
pressure  taps  are  at  off-axis  radii.  The  impeller  centripetal  acceleration  imposes  a  pressure 
gradient  in  the  tube  connecting  the  pressure  tap  to  the  transducer  yielding  a  lower  measured 
pressure  than  the  static  tap  pressure.  A  correction  for  this  difference  is  giver,  in  Equation 
3.1  (Bryan  and  Fleeter,  1987) 


Ptap  =  Paxis  exP  ( )  3  1 

where  fir  and  r  are  the  rotational  speed  and  static  pressure  tap  radius,  and  R  and  T  are  the 
gas  constant  and  temperature,  which  is  assumed  to  be  constant  along  the  tube. 

The  centrifugal  compressor  facility  is  instrumented  with  pressure  probes  for  measuring 
total  and  static  pressure  as  well  as  flow  velocity  and  direction.  Figure  3.2.  The  design  of 
these  probes  is  described  by  Bryer  and  Pankhurst  (1971).  A  five  port  total  pressure  rake 
with  one  cobra  head  is  located  at  the  impeller  inlet  A  three  head  cobra  pobe  is  located  at 
the  impeller  exit  and  two  such  probes  are  located  downstream  of  the  diffuser  vane  row. 
The  inlet  probe  allows  the  total  pressure  and  velocity  to  be  measured  at  five  locations  across 
the  inlet  annulus  as  well  as  the  flow  angle  at  the  mean  radius.  The  three  diffuser  probes 
allow  measurement  of  the  total  pressure,  velocity  and  angle  at  three  positions  across  the 
radial  diffuser. 

The  fluid  temperature  is  measured  at  the  impeller  inlet,  the  diffuser  exit  and  downstream 
of  the  orifice  plate.  This  is  accomplished  using  type  T  thermocouples  in  conjunction  with 
an  HP3497  Type  T  Thermocouple  Interface  Card  with  electronic  ice  point.  The 
thermocouple  voltage  is  measured  using  the  HP3497  voltmeter  and  the  temperature  is 
calculated  by  a  ninth-order  polynomial  specified  by  Omega  Engineering. 

For  the  above  steady  data  acquisition,  the  data  is  sampled  30  times  for  each  pressure 
and  temperature.  The  mean  and  standard  deviation  is  calculated  and  an  uncertainty  is  found 
using  the  t-distribution  with  a  95%  confidence  interval.  The  uncertainties  in  the  pressure 
and  temperature  measurements  are  typically  less  than  1%. 
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3.2  Unstggdy.  Pressure  Instrumentation 

Time  varying  pressures  on  the  diffuser  vanes  are  measured  using  reverse  mounted, 
high  natural  frequency  miniature  pressure  transducers.  Two  diffuser  vanes  are  each 
instrumented  with  eight  Kulite  LQ- 125-5  G  strain  gage  pressure  transducers,  the  pressure 
surface  of  one  vane  and  the  suction  surface  of  the  other.  The  transducers  have  2.18 
mv/kPa  (15  mV/psi)  sensitivity  with  a  70  kHz  natural  frequency.  The  diameter  and 
thickness  are  0.3175  cm  and  0.0889  cm  (0.125  and  0.035  in.)  respectively,  with  a  sensing 
element  diameter  of  0.0965  cm  (0.038  in.).  Two  additional  vanes  are  instrumented  with 
static  pressure  taps  at  the  same  chordwise  locations  as  the  Kulite  transducers  for  time 
averaged  static  pressure  measurement.  These  measurement  locations  are  shown  in  Figure 
3.3. 

The  transducers  are  reverse  mounted  in  a  cavity-pressure  tap  arrangement  that  can  be 
modeled  as  a  Helmholtz  resonator.  The  volume  is  sufficiently  small  to  insure  no  phase 
shift  or  attenuation  within  the  frequency  range  of  interest,  with  accurate  unsteady  pressure 
transmission  at  three  times  the  impeller  blade  pass  frequency.  Tne  signals  from  the  d'ffuser 
vane  mounted  transducers  are  preamplified  using  Analog  Devices’  AD624  precision 
instrumentation  amplifiers  with  a  gain  of  200.  The  signals  are  then  conditioned  using 
adjustable  gain  Ithaco  Model  455  amplifiers.  The  blade  mounted  transducers  are  statically 
calibrated  to  determine  sensitivity  and  to  insure  linearity.  A  dynamic  calibration  (Capece 
and  Fleeter,  1987)  showed  negligible  phase  shift  and  attenuation  to  2100  Hz,  three  times 
the  impeller  blade  pass  frequency. 

The  unsteady  static  pressure  at  the  impeller  exit/diffuser  entrance  is  measured  using 
PCB  Piezotronics  103A  piezoelectric  transducers.  These  transducers  incorporate  self- 
contained  amplifiers  and  have  a  sensitivity  of  217.5  mv/kPa  (1500  mV/psi)  with  a  natural 
frequency  of  13  kHz.  The  diameter  and  thickness  are  0.9525  cm  and  0.5588  cm  (0.375 
and  0.22  in.)  respectively,  with  a  sensing  element  diameter  of  0.254  cm  (0.10  in.).  Two 
such  transducers  are  mounted  through  pressure  taps  to  the  diffuser  tip  wall  at  40.39  cm 
(15.9  in.)  radius  and  33°  angular  spacing. 

The  unsteady  pressure  measurement  signals  are  recorded  using  a  Gould  Biomation 
2805  Waveform  Digitizer.  The  Biomations  are  triggered  by  an  impeller  mounted 
photodetector  and  simultaneously  store  2048  readings  on  four  separate  channels  with  a  user 
adjustable  timebase.  The  raw  buffer  data  is  transferred  to  the  computer  through  a  parallel 
interface.  The  impeller  synchronous  pressure  traces  are  phase  averaged  over  200 
ensembles  to  eliminate  random  unsteadiness  and  isolate  the  periodic  pressure  fluctuations 
which  are  of  primary  interest. 
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3.3  Laser  Doppler  Velocimeter 

The  installation  of  the  Laser  Doppler  Velocimeter  (LDV)  optics  assembly  is  described 
by  Fagan  (1989).  A  dual-beam  LDV  system  operating  in  backscatter  mode  is  utilized, 
schematically  shown  in  Figure  3.4.  In  a  dual-beam  system,  the  beam  emitted  by  the  laser 
is  split  into  two  equal  strength  beams.  The  beams  symmetrically  enter  the  focusing  lens 
and  cross  at  the  lens  focal  plane,  forming  the  probe  volume.  The  light  scattered  from  small 
particles  passing  through  the  probe  volume  is  collimated  by  the  focusing  lens  and  reflected 
into  the  receiving  lens  by  a  mirror.  This  backscattered  light  is  focused  on  the  photodetector 
by  the  receiving  lens. 

A  custom  assembly  has  been  constructed  with  bench  mounted  preliminary  optics. 
Figure  3.5,  connected  via  fiber  optic  link  to  the  final  optics  assembly.  Figures  3.1  and  3.6. 
The  beam  from  a  4  watt  argon-ion  laser  operating  at  514.4  qm  wavelength  is  turned  180° 
by  two  first  surface  mirrors.  A  beamsplitter  divides  the  beams  into  two  equal  power 
beams.  The  beams  are  acousto-optically  frequency  shifted  40  MHz  and  30  MHz 
respectively  using  twin  bragg  cells,  with  an  effective  shift  of  10  MHz.  Frequency  shifting 
is  utilized  to  minimize  fringe  bias  and  eliminate  directional  ambiguity.  The  beams  are  then 
directed  into  a  translator  module  for  proper  positioning  into  the  fiber  optic  input  coupler. 
The  beams  are  transmitted  to  the  final  optics  assembly  through  two  polarization  preserving 
fibers. 

The  final  optics  assembly  contains  the  fiber  optic  output  coupler,  the  focusing  lens,  the 
probe  volume  positioning  minor  and  the  receiving  optics.  The  fiber  optic  output  coupler 
collimates  and  positions  the  beams  leaving  the  fiber  optic  cables.  It  can  be  rotated  360°  by 
a  stepping  motor  for  beam  orientation.  The  beams  emitted  from  the  fiber  optic  output 
coupler  are  directed  into  the  focusing  lens.  The  beams  are  then  reflected  into  the  test 
section  through  a  plexiglass  window  by  the  probe  volume  positioning  mirror,  which  can  be 
rotated  in  two  planes  about  the  beam  axis  by  stepping  motors.  The  LDV  system  operates  in 
backscatter  mode  in  which  the  scattered  light  is  collected  through  the  transmitting  mirror 
and  lens.  The  scattered  light  is  reflected  by  the  positioning  mirror  and  collimated  by  the 
focusing  lens.  It  is  turned  by  a  second  mirror  and  focused  on  the  photodetector  aperture  by 
the  receiving  optics  focusing  lens.  The  entire  final  optics  assembly  can  be  translated  along 
x,  y  and  z  axes  using  stepping  motors,  yielding  a  total  of  six  positioning  axes. 

A  250  mm  focal  length  focusing  lens  is  utilized  in  the  current  investigation,  yielding 
probe  volume  dimensions  of  dm=126.6  pm  (0.005  in.),  hm=127.2  pm  (0.005  in.)  and 
lm=  1266  pm  (0.05  in.),  with  a  volume  of  0.0106  mm3  (6.47  x  10-7  in.3).  The  flow  is 
seeded  using  a  50/50  mixture  of  propylene  glycol  and  ethanol  with  particle  diameters 
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between  0.5  and  1.5  Jim.  The  velocity  ratio  few  these  particles  at  the  impeller  wheel  speed, 
Uw=80m/s,  is  greater  than  98%. 

To  determine  the  probe  volume  position  and  orientation,  an  inverse  ray  tracing 
algorithm  is  used  which  performs  an  exact  ray  trace  from  the  probe  volume,  through  the 
access  window  to  the  focusing  lens.  Prior  to  data  acquisition,  the  six-axis  positioning 
system  is  manually  positioned  on  a  reference  mark  with  proper  orientation.  The  HP1000 
computer  initializes  the  stepper  motor  positioners,  then  the  axis  settings  for  the  desired 
measurement  location  and  direction  are  used  by  the  computer  to  move  the  six-axis 
positioner  to  the  desired  location. 

A  TSI  Model  1990  Counter  Processor  is  utilized  for  photodetector  signal  processing. 
The  Doppler  burst  frequency  is  determined  from  the  time  required  for  a  seed  particle  to 
cross  16  probe  volume  fringes.  Measurements  with  less  than  16  fringe  crossings  are 
invalidated.  Along  with  the  Doppler  burst  frequency  time  measurement,  a  1  MHz  clock 
circuit  is  utilized  to  latch  a  timing  word  at  the  data  ready  signal  from  the  counter  processor 
and  is  reset  at  each  revolution.  This  timing  word  is  transferred  to  the  computer  along  with 
the  LDV  word  and  ured  to  determine  the  impeller  angular  position  at  the  instant  of  the 
measurement.  The  23  bladed  impeller  is  divided  into  460  circumferential  bins,  20  per  blade 
passage,  and  the  data  is  stored  in  the  bin  corresponding  to  the  proper  impeller  position. 
Data  is  collected  at  a  given  location  until  a  requisite  number  of  samples  are  in  each  bin.  In 
this  scheme,  the  time  varying  velocity  signal  around  the  impeller  circumference  can  be 
determined. 

The  particle  velocity  can  be  determined  from  Equation  3.2  (Fagan  and  Fleeter,  1989), 
which  is  an  expression  for  the  photodetector  output  frequency,  0).  The  first  term  on  the 
right  hand  side  of  Equation  3.2  is  the  frequency  difference  of  the  two  illuminating  beams, 
which  have  been  frequency  shifted  in  the  current  investigation.  The  second  term  is 
proportional  to  the  particle  velocity,  u,  and  is  a  function  of  the  beam  half  angle,  y,  which  is 
determined  by  the  ray  tracing  algorithm  using  the  focusing  lens  focal  length,  and  the  beam 
wavelength,  X,  which  is  known. 


.  .  2  sin  y 

CO  =  (0)01  -  0)02)  + - u 

X  3.2 

The  particle  velocity  measurement  is  independent  of  the  photodetector  position. 

Velocity  bias  can  occur  in  LDV  systems  since  a  measurement  can  only  be  made  when  a 
particle  crosses  the  probe  volume.  The  occurrence  of  a  measurement  is  related  to  its 
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instantaneous  velocity,  since  more  particles  will  cross  the  probe  volume  in  a  given  time 
increment  at  greater  velocities.  Roesler,  Stevenson  and  Thompson  (1980)  discovered  that 
constant  frequency  (equal  time)  sampling  at  lower  frequency  than  the  data  realization  rate 
eliminates  the  correlation  between  instantaneous  velocity  and  sampling  rate.  Since  the  data 
in  the  current  investigation  are  stored  in  bins  according  to  impeller  position,  the  data  in  each 
bin  is  effectively  equal  time  sampled  and  velocity  bias  is  eliminated. 


3.4  Compressor  Performance  Analysis 

The  compressor  performance  is  determined  from  the  steady  pressure  and  temperature 
data  and  the  mass  flow  rate  measured  from  the  orifice  plate  static  pressure  taps.  The  total 
and  static  pressures  as  well  as  the  flow  velocities  at  the  impeller  inlet,  impeller  exit,  and 
diffuser  exit  are  determined  from  the  pressure  probe  data.  If  the  density  changes  through 
the  machine  are  small  (<  5%),  which  is  a  good  approximation  in  the  current  investigation, 
the  machine  performance  can  be  expressed  as  the  pressure  rise  coefficient,  VF,  and  the 
efficiency,  q,  as  a  function  of  the  flow  coefficient,  O  (Dixon,  1978). 

The  pressure  rise  coefficients  for  the  impeller  and  the  impeller/diffuser  stage  are  defined 
in  Equation  3.3 


yi=5  PQl..:..?pi 

l/2pU£  3.3a 

vps  =  Pq.2  :  Pqi 

l/2pUw  3.3b 


where  P0i,  Pol.  and  P02  are  the  total  pressures  at  the  impeller  inlet,  impeller  exit  and 
diffuser  exit  respectively,  and  p  and  Uw  are  the  density  and  wheel  speed  at  the  mean 
impeller  exit  radius. 

The  efficiencies  are  defined  as  the  ratio  of  the  useful  energy  imparted  to  the  fluid  to  the 
work  done  by  the  machine,  Equation  3.4. 


Tii  = 


(Pol  ~  P Oi)/p 
UwUei  -  UwiU0i 


3.4a 
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n  _  (Pq2  ~  Poi)/p 

^"UwUei-UwiUei  3.4b 

where  Uei  and  Uei  are  the  absolute  tangential  velocities  at  the  impeller  inlet  and  impeller 
exit  and  Uwi  is  the  wheel  speed  at  the  impeller  inlet  The  flow  coefficient  is  defined  in 
Equation  3.5, 


d>  =  !ia 
Uw 


3.5 


where  Uro  is  the  absolute  radial  velocity  at  the  impeller  exit 

The  pressure  rise  coefficients  and  efficiencies  are  measured  for  the  impeller  and  the 
impeller/diffuser  stage  in  a  vaneless  configuration.  For  vaned  diffuser  configurations,  only 
the  stage  performance  is  determined  due  to  interference  effects  on  the  impeller  exit  probe  by 
the  diffuser  vanes.  The  denominator  of  Equation  3.4  is  determined  for  all  cases  from  the 
vaneless  data  since  the  vaned  diffuser  had  little  effect  on  the  impeller  loading. 


3.5  Impeller  and  Vaned  Diffuser  Steady  Loading  Analysis 
The  static  pressure  distributions  on  the  rotating  impeller  blades  and  the  diffuser  vanes 
were  measured.  A  pressure  and  pressure  difference  coefficient  for  the  impeller  are  defined 
in  Equations  3.6  and  3.7. 


Cp  = 


P pres  ~  P suet 

l/2pu£ 


3.6 

3.7 


This  steady  pressure  coefficients  are  presented  as  a  function  of  the  nondimensional  radius 
along  the  blade,  r/r0,  where  r0  is  the  impeller  exit  radius. 

The  pressure  coefficient  and  the  pressure  difference  coefficient  for  the  diffuser  vanes  is 
given  in  Equations  3.8  and  3.9. 


CP  = 


— E_ 
l/2pUj 


3.8 


where  Uj  is  the  steady  absolute  velocity  at  die  diffuser  entrance. 

The  pressure  coefficients  are  presented  as  a  function  of  the  nondimensional  distance 
along  the  blade  chord,  x/c.  The  steady  lift  coefficient  is 


3.10 


where  c  is  the  vane  chord. 

Hie  diffuser  vane  mean  incidence  angle  is,  I  =  o  -  <*i,  where  o  is  the  diffuser  vane  row 
stagger  and  ai  is  the  vane  leading  edge  absolute  flow  angle. 

The  uncertainties  can  be  calculated  for  die  performance  and  loading  parameters  from  the 
uncertainties  of  the  raw  data.  Typically,  die  uncertainties  in  the  above  coefficients  are  less 
than  2%. 


3.6  LPV  Velocity  Analysis 

The  LDV  data  arc  sampled  at  each  measurement  location  until  approximately  2000 
samples  are  taken  for  each  blade  passage.  As  previously  mentioned,  the  data  are  stored  in 
460  bins  spaced  circumferentially  around  the  impeller,  20  per  blade  passage.  The  mean 
and  standard  deviation  are  calculated  for  the  data  in  each  bin,  and  the  total  unsteady  velocity 
distribution  around  the  impeller  is  constructed  from  the  mean  velocity  in  each  bin.  Fagan 
(1989)  has  shown  that  the  uncertainty  of  the  mean  in  each  bin  is  approximately  0.15%. 

The  time  averaged  radial  and  tangential  velocity  components,  Ur  and  Ug,  are  evaluated 
by  taking  the  average  of  the  radial  and  tangential  velocity.  A  Fast  Fourier  Transform  (FFT) 
is  performed  on  the  velocity  data  to  determine  the  harmonic  content  of  the  unsteady  portion 
of  the  data.  The  first  harmonic  unsteady  radial  and  tangential  velocities,  vr  and  vg,  are 
determined  from  the  FFT. 
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3.1.  Unsteady-Pressure  Data  Analysis 

The  impeller  synchronous  ensemble  averaged  unsteady  airfoil  static  pressure  and 
impeller  exit  static  pressure  are  measured  using  the  high-response  transducers  mentioned 
previously.  The  ensemble  average  eliminates  random  unsteadiness  and  preserves  the 
pressure  variations  which  are  synchronous  with  the  impeller  rotation.  The  measured 
voltage  signals  are  converted  to  pressure  through  the  transducer  calibration  and  an  FFT  is 
performed  on  the  pressure  signals  to  determine  the  first  harmonic  amplitude  and  phase. 
The  first  harmonic  unsteady  pressures  can  then  be  expressed  in  complex  form, 
p  =  p  exp(-icot) ,  where  p  is  the  complex  amplitude  of  the  first  harmonic  unsteady  pressure. 

An  airfoil  surface  unsteady  pressure  coefficient  is  defined  by 


Pvm  Ui  3.11 

where  Ui  is  the  steady  absolute  velocity  at  the  diffuser  entrance  and  v„i  is  the  transverse 
velocity  which  is  the  unsteady  first  harmonic  velocity  at  the  airfoil  leading  edge  normal  to 
the  mean  flow 


vnl  =  vrjsin  <Xi  -  veicos  ai 


3.12 


where  ai  is  steady  absolute  flow  angle  at  the  diffuser  entrance. 
A  pressure  difference  coefficient  is  defined  as 


CAP  = 


PU-PL 

pvm  U! 


where  pu  and  pL  are  the  upper  and  lower  surface  unsteady  pressures. 
The  complex  unsteady  lift  coefficient  is  given  in  Equation  3.14 


CL  = 


dx/c 


3.13 


3.14 


This  lift  coefficient  expresses  the  unsteady  lift  based  on  the  leading  edge  steady  and 
unsteady  normal  velocities.  A  second  lift  coefficient  based  on  the  unsteady  pressure 
amplitude  and  steady  velocity  is  given  in  Equation  3.15 
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CU» 


JEUiEiW 

l/2pUw 


3.15 


This  coefficient  indicates  the  vane  loading  level  as  a  function  of  the  machine  operating 
point 

The  reduced  frequency  and  interblade  phase  angle  are  defined  as 


k 


Ui 


3.16 


0b”-2it 


Nt 

Nv 


3.17 


where  Qb  is  the  impeller  blade  pass  frequency,  and  Nb  and  Nv  are  the  number  of  impeller 
blades  and  number  of  diffuser  vanes  respectively. 


3,8,  Summary 

The  Purdue  Research  Centrifugal  Compressor  instrumentation  allows  extensive  steady 
and  time-varying  measurements  to  be  made  in  the  inlet  impeller  and  diffuser.  The  steady 
static  and  total  pressure  along  with  flow  direction  and  velocity  at  the  impeller  inlet  exit  and 
diffuser  exit  are  measured  along  with  static  pressures  on  the  channel  walls,  diffuser  vanes 
and  impeller  blades.  The  time  varying  static  pressures  are  measured  on  the  diffuser  vanes 
and  diffuser  endwalls  at  the  impeller  exit  The  detailed  two-dimensional  time-varying 
diffuser  velocity  field  is  measured,  with  the  capability  of  three-dimensional  measurements 
throughout  the  axial  inlet  impeller  podges  and  vaned  diffuser. 


Figure  3.1  Rotating  Scanivalve  and  LDV  Final  Optics 


Vane  Leading  Edge  1.5 


Figure  3.3  Diffuser  Vane  Unsteady  Pressure  Measurement  Locations 


Photodetector 
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Figure  3.4  Dual  B-  am  LDV  System 
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Figure  3.5  LDV  Table-Mounted  Preliminary  Optics 
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Figure  3.6  LDV  Final  Optics  Assembly 
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CHAPTER  4 

RADIAL  DIFFUSER  WAKE  ANALYSIS 

Several  analyses  have  been  developed  to  predict  the  unsteady  pressure  distribution  on 
axial  flow  cascades  (Whitehead,  1960;  Whitehead  1978;  Chiang  and  Fleeter,  1988)  and 
recently  in  radial  flow  cascades  (Bryan  and  Fleeter,  1990).  Typically,  these  unsteady 
aerodynamic  models  are  linearized,  i.e.,  the  unsteady  velocities  and  pressures  are 
considered  small  perturbations  superimposed  on  a  steady  mean  flow,  with  second  order 
and  higher  terms  neglected.  Hie  unsteady  velocity  field  is  predicted  by  solving  for  the 
unsteady  velocity  potential,  with  boundary  conditions  given  far  from  the  blade  row  and 
also  on  the  blade  surface.  When  solving  for  the  unsteady  pressure  distribution  generated 
by  an  upstream  wake,  the  velocity  induced  by  the  wake  on  the  airfoil  surface  is  specified  as 
a  boundary  condition.  In  most  models,  the  wake  is  assumed  to  be  convected  through  the 
cascade,  undisturbed  by  the  airfoils.  Only  the  effect  of  the  wake  on  the  airfoil  unsteady 
pressure  distribution  is  considered. 

In  this  chapter,  a  mathematical  model  is  developed  to  determine  the  behavior  of  vortical 
wakes  in  a  radial  vaneless  diffuser  for  use  in  unsteady  aerodynamic  models  for  radial  flow 
blade  rows.  Several  investigations  have  considered  both  steady  and  time-variant  radial 
diffuser  flows,  including  unsteady  wake  behavior  (Dean  and  Senoo,  1960;  Senoo  and 
Ishida,  1975).  These  unsteady  wake  models  are  control  volume  solutions  and  are  not 
directly  applicable  for  use  as  boundary  conditions  in  unsteady  airfoil  models.  Therefore,  a 
model  is  developed  which  describes  the  behavior  of  small  amplitude  vortical  impeller 
wakes  convected  through  the  radial  diffuser  by  the  mean  flow. 

4J.Analy.sis 

The  flowfield  in  a  parallel  wall  radial,  vaneless  diffuser  is  considered.  The  flow  in  the 
radial  diffuser  is  assumed  to  be  incompressible  and  inviscid.  The  imposed  boundary 
conditions  are  derived  from  the  flow  behind  a  radial  impeller,  with  these  boundary 
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conditions  assumed  to  be  equal  to  the  radial  diffuser  inlet  boundary  conditions  where 
velocity  measurements  are  made. 

The  radial  and  tangential  velocities  are  considered  to  consist  of  a  small  amplitude 
unsteady  velocity  superimposed  on  the  mean  steady  velocity. 

vr  =  Ur  +  vr  4.1a 

ve  =  Ue  +  v0  4.1b 

The  nonuniform  circumferential  flow  leaving  the  relative  reference  frame  of  the  centrifugal 
impeller  is  unsteady  in  the  absolute  reference  frame  of  the  radial  diffuser.  Figure  4.1. 
Since  the  flow  is  inviscid,  the  flowfield  can  be  described  using  the  unsteady  Euler 
equations.  Assuming  the  unsteady  flow  to  be  a  small  disturbance  superimposed  on  the 
steady  mean  flow,  Equations  4.1,  the  steady  and  unsteady  velocities  can  be  described  by 
separate  equations,  where  second  order  terms  in  the  unsteady  velocities  are  neglected. 


(  V-V)  Qs  =  0 

4.2a 

v  V»  =  -ft* 

4.2b 

a  -  \ 

-+  V-V  Q  =  0 

.a*  / 

4.3a 

< 

ii 

i 

s 

4.3b 

where  V  is  the  steady  velocity  vector,  Q  and  Qj  are  the  unsteady  and  steady  vorticity  and  \jr 
and  \jrs  are  the  unsteady  and  steady  stream  functions.  The  steady  flow,  Equations  4.2,  is 
independent  of  the  unsteady  flow,  however,  the  unsteady  flow,  Equations  4.3,  is 
dependent  on  the  steady  flow. 

The  boundary  conditions  for  Equations  4.2  and  4.3  in  the  relative  reference  frame  of 
the  impeller  are  depicted  in  Figure  4.1  and  given  by  Equations  4.4. 


vr  =  W0  cos  Po  +  vro  exp  (i  Nb  0r)  @  r=r0 


4.4a 
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v0  =  W0  sin  p0  +  exp  (i  Nb  0r)  @  r=r0  4.4b 

vr  =  v0  =  0  as  r-)oo  4  4C 

where  W0  is  the  mean  velocity  in  rotating  frame,  Nb  is  the  number  of  impeller  blades,  po  is 
the  mean  relative  flow  angle,  r0  is  the  impeller  exit  radius  and  vro  and  V0o  are  the  complex 
amplitudes  of  the  wake  velocities  at  the  impeller  exit 

The  boundary  conditions  specified  in  Equations  4.4  can  be  transformed  to  the 
stationary  reference  frame  by  adding  the  wheel  speed  Uw  to  the  tangential  velocity  and 
expressing  the  relative  tangential  coordinate,  0r,  in  terms  of  the  absolute  tangential 
coordinate  0,  and  time,  t  The  absolute  coordinate  0,  is  related  to  the  relative  coordinate  0r, 
by  0  =  0r  +  tot .  Substituting  this  into  Equations  4.4  and  adding  the  wheel  speed  yields  the 
radial  and  tangential  velocity  boundary  conditions  at  the  impeller  exit 

vr  =  W0  cos  Po  +  Vjo  exp  (ik(0  - 1 ))  @  r  =  r0  4.5a 

v0  =  W0  sin  Po  +  exp  (ik(0- 1 ))  +  Uw  @  r  =  r0  4.5b 

where  0  =  0/(U/Uw)  and  t  is  scaled  by  Tq/Uo,  where  U0  is  the  steady  absolute  velocity  at 
the  impeller  exit 

The  solution  to  Equations  4.2,  assuming  no  variation  in  the  steady  velocity  with  0,  are 

Ur~  Urn/  r  4.6a 

Ue  =  U0o/r  4.6b 

where  Uro  and  Ueo  are  the  steady  mean  absolute  velocities  at  the  impeller  exit  radius,  r0, 
and  r  is  the  nondimensional  radius,  r/r0. 

The  unsteady  vorticity  equation,  Equation  4.2a,  in  radial  coordina.es  is 

i+Utitiki_)n=o 
at  ar  r  a©/ 


4.7 
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With  Ur  and  Uq  determined.  Equation  4.7  can  be  solved  for  the  unsteady  vorticity,  £2.  The 
unsteady  vorticity  is  considered  an  harmonic  function  of  time, 


£2  =  £2  exp  (-i  to  t) 


4.8 


where  £2  is  the  complex  unsteady  vorticity  amplitude  which  is  a  function  of  r  and  0. 

Substituting  Equations  4.4  and  4.8  into  Equation  4.7  yields  the  following  equation  for 
the  unsteady  vordcity  amplitude,  Q. 


r  dr  r 2  30/ 


£2  =  0 


4.9 


where  Uro  and  Ueo  are  nondimensionalized  by  the  absolute  velocity  leaving  the  impeller, 
U0,  and  k  is  the  reduced  frequency,  coro/Uo 

Equation  4.9  is  a  first  order,  linear  partial  differential  equation  with  the  following 
solution  (Zachmanoglou  and  Thoe,  1986) 


£2  =  £2o  expMk-  r*)  fir  exp(,  ^  0))  4. 10 

l  Uro  Uqq 


where  £2q  is  the  complex  amplitude  of  the  unsteady  vorticity  and  f  is  a  functional  form  to  be 
determined. 

The  functional  form,  f,  can  be  determined  from  the  boundary  conditions.  Equations 
4.5.  The  vorticity  is  defined  as 


^  _  l  3(rve)  idvr 

r  *  "  r39 


4.11 


Since  £2  is  determined  from  partial  derivatives  of  the  velocity,  the  dependence  of  £2  on  0 
must  be  the  same  as  the  dependence  of  the  velocity  derivatives,  which  is  an  exponential. 
Therefore,  the  functional  form,  f,  must  be 

fix)  =  x  ^NbfUeo/Uro))  412 


The  unsteady  vorticity,  determined  from  Equation  4.8,  becomes, 
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Q  =  Qoexp(A  In  r  +  B  r2)  exp(ik6 ) 


where 


A  =  -ik^a^Q-  =  -ik - -  and  B=Jk_  =  — Ul_ 

UwUro  (tan  Oq  -  tan  Po)  cos  Oq  2Uro  2cosoto 


and  ao  is  the  steady  mean  flow  angle  in  the  absolute  coordinate  system. 

To  determine  the  unsteady  stream  function,  V,  the  expression  for  the  unsteady 
vorticity.  Equation  4.13,  is  substituted  into  Equation  4.3b. 


— +  =  -ro2Oo«p(A  In  r  +  B  r2)  exp(ik0) 

dr2  r  dr  30  2 


The  solution  to  Equation  4.14  is  determined  as  a  complex  Fourier  series  by  a  finite  Fourier 
transform.  Equation  4.15,  taking  advantage  of  the  fact  that  the  solution  is  periodic  in  ft 


VM)=  E  c„(r)e*« 


The  left  hand  side  coefficients  are  defined  by 


-tf 


y(r,0)  e-‘n*d0 


and  the  source  term  coefficients  are 


Q(r,0)e’in9d0 
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Substituting  Equations  4.16  and  4.17  into  Equation  4.14  yields 


(£F  C  =  V.%exp(Ator  +  Br2> 


Equation  4.18  is  a  Cauchy-Euler  differential  equation.  The  complementary  solution  is 
given  in  Equation  4.19  (Ross,  1980). 

Cn  =  ai  rk2  +  a2  r‘k2  4'  ^ 

where  k2  =  k(Uo/Uw). 

The  particular  solution  to  Equation  4.19  is  found  by  variation  of  parameters.  The 
complete  solution  for  cn  is  given  in  Equations  4.20  and  4.21. 

c„  =  ai/k2  +  a2  rk 2  +  r<^Qp-rk2  j  exp(A  r  +  Br2)  r  1  -  k2  dr 


.  £s£-Qa.  r-k2  f 

2  k2  l 


exp(A  In  r  +  B  r2)  r1  +  k2  dr 


c.n  =  bj  rk2  +  b2  r_k 2 


The  stream  function  can  be  constructed  from  the  Fourier  series.  Equation  4. 15. 


v  =  rQ2  pQ.  (rk2  ( 

M  J, 


exp(A  In  r  +  Br2)  r 1  *  k2  dr  - 


r*k2  I  exp(A  In  r  +  B  r2)  r 1  +  k2  dr  exp(ik(6  - 1 ))  + 


(ai  rkz  +  a2  r*2)  exp(ik(0  - 1 ))  +  (bi  rkz  +  b2  r^*)  exp(ik(-0  - 1 )) 
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Since  the  solution  is  not  a  function  of  e**®,  bi  and  b2  are  zero.  The  unsteady  radial  and 
tangential  velocities  are  evaluated  from  the  stream  function. 


vr 


r39  2 


exp(A  In  r  +  Br2)  r1  *  k2  dr  - 


r-k2- 1 


exp(A  In  r  +  B  r2)  r1  +  k2  dr 

i 


exp(ik(0  -t))  + 


(i^_  rk2  - 1  +  i^2.  rk2  -  l)exp(ik(0  - 1 )) 


and 


v0 


=  .^  =  .£a|k  Lk2-1  I  exp(A In r  +  Br2) r 1  •  k2  dr 

2  \  J\ 

-■f 


rk2  - 1  I  exp(A  In  r  +  B  r2)  r1  +  k2  dr 


exp(ik(0- 1))  + 


(' ^ rk2'l+  ^ rkl ' exp(ik(6 - 1 )) 


4.24 


The  unsteady  wake  velocities,  Equations  4.23  and  4.24,  are  composed  of  a  vortical 
velocity  component  and  an  iirotational  velocity  component.  The  irrrotational  component 
corresponds  to  the  homogeneous  solution  of  Equation  4.14,  with  the  rotational  component 
corresponding  to  the  particular  solution,  which  contains  the  Qq  term  in  Equations  4.23  and 
4.24. 

The  unsteady  wake  velocity  solution,  Equations  4.23  and  4.24,  is  determined  by 
specifying  Oq,  p0  and  k.  The  integrals  are  evaluated  numerically  using  Simpson’s  1/3 
rule.  The  constants  at  and  a2  are  evaluated  from  the  boundary  conditions,  Equations  4.3a 
and  4.5b.  The  vorticity  magnitude,  fio,  is  evaluated  from  the  third  boundary  condition, 
Equation  4.4c. 

The  unsteady  wake  velocities,  vr  and  vq,  are  for  a  vaneless,  parallel  wall,  radial 
diffuser.  The  steady  flowfxeld  is  determined  by  the  absolute  flow  angle,  cXq.  The  wake 
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propagation  is  determined  by  the  relative  flow  angle,  p0,  the  reduced  frequency,  k,  and  the 
boundary  conditions  of  the  wake,  vro  and  ve0 .  In  principle,  the  wave  number,  k2,  can 
assume  any  value  for  given  Oq,  po  and  k.  In  practice,  k2  is  equal  to  the  number  of  impeller 
blades,  Nb,  and  will  always  be  an  integer.  Thus,  given  the  number  of  impeller  blades,  Nb, 
the  relative  flow  angle,  Po,  and  the  absolute  flow  angle,  Oq,  the  unsteady  wake  velocities 
can  be  determined  from  Equations  4.23  and  4.24. 


4.2-Rcsults 

The  unsteady  wake  model  developed  herein  is  utilized  to  demonstrate  the  effect  of 
steady  centrifugal  compressor  parameters,  including  the  absolute  flow  angle,  a0,  the 
relative  flow  angle,  P0,  and  the  number  of  impeller  blades,  Nb.  The  unsteady  wake 
velocities,  vr  and  vg,  are  constant  amplitude  harmonic  functions  of  the  tangential  coordinate 
9,  *>nd  time  r ,  at  a  particular  value  of  r.  However,  the  velocity  varies  in  amplitude  with  r 
due  to  the  radial  geometry  and  thus  variation  in  steady  velocity  with  radius.  In  the 
following  figures,  the  radial  dependence  of  vr  1  vg  is  presented  as  a  function  of  the 
nondimensional  radius,  assuming  no  variation  of  the  relative  flow  angle,  f30,  with  the 
tangential  direction,  0.  A  parametric  study  is  performed  by  varying  a o,  Po  and  Nb-  The 
baseline  values  are  oto  =  60°,  Nb  =  12  and  P0  =  -45°.  Negative  p0  values  indicate  impeller 
blades  leaning  in  the  opposite  direction  of  rotation. 

Figures  4.2  and  4.3  show  the  results  for  the  baseline  case  with  aQ  =  60°  and  for 
a o  =  45°.  The  radial  velocity  monotonically  approaches  zero  whereas  the  circumferential 
velocity  increases  to  a  maximum  before  approaching  zero.  This  behavior  is  typical  for 
Po  <  -45°.  The  radial  velocity  decreases  more  rapidly  than  the  tangential  velocity,  with 
both  components  decreasing  more  rapidly  as  the  absolute  flow  angle,  Oq  is  increased. 

Figures  4.4  and  4.5  show  the  wake  velocities  for  p0  =  -30°  and  p0  =  0°.  Changing 
Po,  the  relative  flow  angle,  changes  the  boundary  conditions  for  vro  and  vg0>  under  the 

assumption  of  circumferentially  constant  Pq.  For  various  values  of  P0,  the  radial  velocity 
component  demonstrates  behavior  similar  to  the  baseline,  with  only  a  magnitude  variation. 
The  tangential  velocity  again  increases  before  its  eventual  decay,  but  the  maximum  occurs 
at  larger  values  of  r  for  Po  closer  to  zero. 

Figure  4.6  presents  the  wake  velocities  with  Nb  =  20.  Variation  in  the  number  of 
impeller  blades,  Nb,  has  a  minimal  affect  on  the  amplitude  of  the  wake  velocities,  with  only 
the  phase  changing  more  rapidly  with  increasing  Nb. 
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The  predicted  wake  velocities,  vr  and  vq,  will  be  compared  to  the  measured  wake 
velocities  in  the  vaneless  radial  diffuser  configuration  of  the  Purdue  Research  Centrifugal 
Compressor.  As  mentioned  previously,  these  velocity  expressions  will  also  be  utilized  as  a 
boundary  condition  in  the  unsteady  radial  airfoil  cascade  analysis  in  Chapter  5. 


CHAPTER  5 

RADIAL  CASCADE  UNSTEADY  AIRFOIL  THEORY 


Much  attention  has  been  devoted  in  the  literature  to  thr  unsteady  aerodynamics  of 
isolated  airfoils  and  axial  flow  cascades.  In  contrast,  relatively  little  attention  has  been 
directed  at  unsteady  flow  in  radial  cascades,  with  the  experimental  investigations  of  Inoue 
and  Cumpsty  (1984)  and  Arndt  et  al.  (1988)  two  of  the  few  examples.  Radial 
turbomachine  cascade  flow  Helds  and  geometries  are  different  than  those  for  axial  flow 
machines.  Due  to  these  differences,  axial  flow  cascade  unsteady  aerodynamic  models  are 
not  directly  applicable  to  radial  flow  cascades.  Thus,  unsteady  aerodynamic  models 
applicable  to  radial  flow  airfoil  cascades  are  needed. 

An  analysis  is  developed  to  predict  the  unsteady  fluid  dynamics  of  radial  flow  cascades 
while  taking  advantage  of  the  existing  research  for  axial  flow  turbomachines.  This  is 
accomplished  by  developing  a  conformal  transformation  which  allows  a  potential  flow 
solution  for  a  given  axial  flow  cascade  to  be  applied  to  a  radial  cascade  (Smith,  1970; 
Faulders,  1956).  Thus,  this  transformation  mathematically  unwraps  the  steady  and 
unsteady  flowfields  in  a  radial  cascade  into  a  corresponding  axial  cascade  flow.  This 
transformation  is  then  coupled  to  an  incompressible  flat  plate  cascade  analysis  (Whitehead, 
1960;  1978).  The  effects  of  the  cascade  geometry  on  the  unsteady  fluid  dynamics  of  radial 
cascades  are  then  demonstrated  by  applying  this  model  to  a  series  of  radial  flow  cascade 
geometries. 


5.1  Transformation  Analysis 

For  a  two-dimensional,  incompressible  and  irrotational  flow,  the  velocity  potential  <|>  is 
described  by  the  Laplace  equation. 
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From  the  theory  of  complex  variables,  any  function  <|>  of  the  cartesian  variables  x  and  y 
which  satisfies  the  Laplace  equation  can  be  expressed  solely  in  terms  of  the  complex 
variable  z  =  x  +  iy.  Also,  a  conjugate  function  y  exists  such  that  the  complex  potential, 
w  =  <}>  +  iy  is  a  function  only  of  the  complex  variable  z,  with  $  and  y  individually 
satisfying  the  Laplace  equation  (Churchill,  1984). 

The  converse  is  also  true.  Given  any  function  of  a  complex  variable  £  =  r\  +  i£  =  f(z), 
the  real  and  imaginary  components  of  £  individually  satisfy  the  Laplace  equation  provided 
£  is  an  analytic  function  of  z.  Any  complex  function  of  a  complex  variable, 
C  =  T](x,y)  +  i£(x,y)  =  f(x+iy)  can  also  be  considered  a  transformation  from  the  variables  x 
and  y  to  the  variables  TJ  and  In  this  case,  the  transformation  is  termed  a  conformal 
transformation. 

If  a  solution  to  the  Laplace  equation  can  be  found  for  a  given  geometry,  the  solution  in 
any  other  geometry  can  be  found  by  applying  an  appropriate  conformal  transformation. 
For  a  radial  cascade,  the  simple  transformation  given  in  Equation  5.2  transforms  the 
solution  to  a  two-dimensional,  infinite,  axial  flow  cascade,  as  depicted  in  Figure  5.1. 

£  =  lnz  5.2 

where  z  is  the  coordinate  in  the  radial  diffuser  and  £  is  the  coordinate  in  the  transformed 
plane.  Thus,  this  transformation  unwraps  the  steady  and  unsteady  radial  cascade  flow 
fields  into  those  of  a  corresponding  axial  cascade.  This  allows  a  potential  flow  solution  for 
a  given  axial  cascade  to  be  applied  to  a  corresponding  radial  cascade. 

Steady  potential  flow  streamlines  in  a  parallel  wall  radial  diffuser  follow  logarithmic 
spirals  as  shown  in  Chapter  1  and  given  in  Equation  5.3. 

In  Gr-)  =  — B — 

r2  tan  aL  5.3 

where  aL  is  the  angle  between  the  streamline  and  a  radial  line,  i.e.,  ol  is  the  logarithmic 
spiral  angle,  r  denotes  the  radius,  T2  is  a  reference  radius  and  6  is  the  angular  coordinate. 
These  streamlines  are  transformed  into  parallel  straight  streamlines  by  Equation  5.2.  Thus, 
the  flat  plate  airfoils  in  axial  flow  unsteady  aerodynamic  models  are  transformed  into 
logarithmic  spiral  airfoils  in  the  radial  cascade. 
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A  modified  form  of  the  transformation  specified  by  Equation  5.2  is  used  for  the  actual 
conformal  transformation.  Namely,  the  transformation  given  in  Equation  5.4  is  utilized  to 
maintain  the  Kutta  condition  through  the  transformation  in  the  unsteady  flow. 

£  =  r2(lnz-lnr2)  54 

where  z  is  the  coordinate  for  the  radial  geometry  and  £  is  the  coordinate  for  the  transformed 
axial  geometry.  Equation  5.4  expressed  in  radial  coordinates  is, 

n  +  i£  =  f2(ln  (£)  +  i0)  5  5 


The  velocities  in  the  two  coordinate  systems  are  related  by  Equation  5.6. 


Un-iU$ 

hi 


=  Ur  -  iUe 


5.6 


where  and  U§  are  the  axial  cascade  velocity  components  in  the  axial  and  tangential 
directions  respectively,  with  Ur  and  Ue  denoting  the  radial  cascade  velocity  components  in 
the  radial  and  tangential  directions. 

The  Kutta  condition,  applied  to  both  the  steady  and  unsteady  flow  fields,  is  maintained 
through  the  transformation.  In  a  steady  incompressible,  inviscid  flow,  the  pressure  at  any 
point  in  the  flowfield  is  given  by  the  following  Bernoulli  equation. 


E+u!= 

P  2 


5.7 


Since  the  right  hand  side  of  Equation  5.7  is  constant,  the  steady  pressure  is 
proportional  to  the  square  of  the  steady  velocity.  Kutta’s  hypothesis  requires  that  the 
pressure  be  continuous  at  the  airfoil  trailing  edge  and  that  the  trailing  edge  velocity  is  finite. 
Therefore  the  velocities  on  the  upper  and  lower  airfoil  surfaces  must  approach  a  common 
value  at  the  trailing  edge.  Since  the  upper  and  lower  velocities  are  scaled  by  the  same  xlxi 
as  appears  in  the  radial  coordinate  form  of  the  transformation.  Equation  5.5,  the  Kutta 
condition  requirement  of  no  steady  pressure  discontinuity  at  the  airfoil  trailing  edge  is 
satisfied  automatically  independent  of  the  value  of  tj. 
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The  unsteady  pressure  is  described  by  the  following  incompressible  unsteady  Bernoulli 
equation. 


P  2  St 


5.8 


In  axial  flow  cascade  linearized  aerodynamic  models,  the  unsteady  flow  is  considered 
to  be  an  harmonic  small  perturbation  to  a  steady  uniform  mean  flow  in  the  chordwise 
direction.  In  the  chordwise  direction,  the  pressure,  velocity  components,  and  velocity 
potential  are  thus  given  by  Equation  5.9. 


p  =  P  +  p  exp  (-i  0)  t) 

5.9a 

u  =  U  +  u  exp  (-i  oo  t) 

5.9b 

v  =  v  exp  (-i  oo  t) 

5.9c 

<)>  =  4>  +  <j>  exp  (-i  oo  t) 

5.9d 

where  P,  <D,  and  U  are  the  steady  values  of  the  pressure,  velocity  potential  and  freestream 
chordwise  velocity,  with  p,  <|>,  u  and  v  the  complex  amplitudes  of  the  unsteady  pressure, 
velocity  potential,  and  chordwise  and  normal  velocity  components  which  are  functions  only 
of  the  spatial  coordinates.  Note  that  in  the  axial  flow  cascade,  U  and  P  are  constant.  In  the 
radial  cascade,  U  is  proportional  to  1/r  and  P  is  proportional  to  1/r2. 

The  linearized  unsteady  Bernoulli  equation.  Equation  5.10,  is  determined  by 
substituting  Equation  5.9  into  Equation  5.8  and  retaining  only  first  order  terms  in  the 
perturbation  quantities. 


—  +  uU  -  iocxJ>  =  F(t) 

P  5.10 

Applying  die  ’linearized  unsteady  Bernoulli  equation  to  the  upper  and  lower  surfaces  of 
a  reference  airfoil  in  either  the  axial  or  radial  cascade  results  in  the  following  equation  for 
the  complex  unsteady  pressure  difference  across  the  reference  airfoil  chordline. 

^B.^.uCuu-uO  +  iaX^u-W 


5.11 


The  value  of  the  velocity  potential  0  at  each  point  in  the  axial  cascade  is  the  same  as  at 
its  image  point  in  the  corresponding  radial  cascade.  However,  the  axial  cascade  velocity 
and  velocity  components  must  be  multiplied  by  the  metric  of  the  transformation,  l/(r/r2). 
Thus,  the  first  term  on  the  right  hand  side  of  Equation  5.1 1  must  be  multiplied  by  l/(r/r2)2. 
For  the  pressure  to  be  continuous  and  therefore  the  pressure  difference  across  the  airfoil  to 
vanish  at  the  trailing  edge,  rfa  must  be  equal  to  unity  at  the  trailing  edge,  i.e.,  the  radial 
cascade  airfoil  trailing  edge  is  defined  by  r  =  t^. 


5.2  Cascade  Model 

For  a  flat  plate  axial  flow  cascade  at  zero  incidence  in  an  incompressible  flow,  the 
parameters  necessary  to  specify  the  problem  are  the  solidity,  c/s,  the  reduced  frequency, 
k  =  coc/U,  the  stagger  angle,  a,  and  the  interblade  phase  angle,  0b-  For  the  radial  cascade, 
the  important  parameters  are  the  number  of  radial  airfoils,  Nv,  the  airfoil  radial  extent 
divided  by  the  trailing  edge  radius,  Ar/r2,  the  reduced  frequency  based  on  the  airfoil  leading 
edge  velocity,  k  =  coc/U  i,  the  stagger  angle,  O,  which  is  equal  to  the  logarithmic  spiral 
angle,  <Xl,  and  the  interblade  phase  angle,  0b- 

Figure  5.2  shows  the  flat  plate  axial  flow  cascade  transformation  of  a  radial  cascade. 
To  apply  existing  axial  flow  cascade  unsteady  aerodynamic  analyses  to  radial  cascades,  it  is 
necessary  to  determine  the  axial  cascade  parameters  which  correspond  to  a  specified  radial 
cascade  geometry.  Given  the  parameters  defining  the  radial  cascade,  the  appropriate 
corresponding  parameters  for  the  axial  cascade  unsteady  aerodynamic  models  are  calculated 
through  the  transformation  specified  in  Equation  5.5.  In  the  following  equations, 
unsubscripted  values  refer  to  the  radial  cascade  and  values  with  the  A  subscript  refer  to  the 
axial  cascade  transformation  of  the  radial  cascade. 

The  chord  in  the  transformed  axial  cascade  is  given  by  Equation  5. 12 

_  -Ti  In  (1  -  Ar/r2) 

A  coso  5.12 


The  transformed  axial  cascade  solidity,  ca/sa  is 
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The  relationship  between  the  reduced  frequencies  for  the  axial  and  radial  cascades  is  given 
by 


kA  =  _k  In  (1  -  Ar/T2> 

(Ar/r2)(l  -  Ar/r2)  5.14 

where  k  is  the  radial  cascade  reduced  frequency. 

The  state-of-the-art  axial  flow  cascade  unsteady  aerodynamic  models  predict  the 
unsteady  aerodynamic  lift  and  moment  as  well  as  the  pressure  distribution  for  airfoils 
executing  translational  and  torsional  motions  and  subjected  to  convected  harmonic  vortical 
gusts.  The  unsteady  pressure  difference,  and  unsteady  lift  and  moment  coefficients  for 
translation,  torsion  and  gust  are  defined  as 

r  AP  Ap  Ap 

'-'AP  .  »  -  » 

PU h  pU 2a  pUv  5.15a 

CL  =  -Lift- t  -Lift  -,  Lift.. 

pU/ic  pU2ac  pUvc  5.15b 

CM  =  Moment ,  Moment  t  Moment 

pUAc2  pU2ac2  pUvc2  5 15c 

where  h  and  a  are  the  translation  and  torsion  airfoil  motions  h-ho  exp  (-i  k  t) , 
a=Oo  exp  (-i  k  t)  and  v  is  the  wake  or  gust  velocity  at  the  airfoil  leading  edge, 
v  =  v„  exp  (-i  k  t) .  In  the  previous  expressions,  ho,  Oq  and  v0  are  the  motion  and  velocity 
amplitudes.  In  the  radial  cascade,  U  and  v  are  the  values  at  the  airfoil  leading  edge. 

The  unsteady  pressure  distribution  for  the  radial  cascade  is  calculated  from  the  unsteady 
velocities  in  the  axial  cascade  through  the  transformation  defined  by  Equations  5.5  and  5.6, 
with  the  unsteady  aerodynamic  lift  and  moment  coefficients  calculated  from  this  by 
integration. 


5.3  Boundary  Condition  Transformation 

The  radial  cascade  logarithmic  spiral  airfoils  undergoing  translation  and  torsion  motions 
do  not  transform  exactly  to  flat  plates  in  the  axial  cascade  undergoing  the  same  motions. 
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To  correct  for  this,  the  airfoils  in  the  radial  cascade  are  assumed  to  undergo  small  amplitude 
torsional  and  translation  motions,  Figure  5.3,  with  the  corresponding  motions  in  the  axial 
cascade  calculated  through  the  transformation,  Equation  5.5.  The  translation  and  torsion 
boundary  conditions  in  the  form  of  upwash  velocities  for  the  radial  cascade  are  given  in 
Equations  5.16  and  5.17. 


5.16 


va  =  Ua  +  ijjZ(x-xta) 


5.17 


where  x  is  the  distance  along  the  vane  nondimensionalized  by  the  chord  and  Xea  is  the 
nondimensional  elastic  axis  location  defined  in  Figure  5.2. 

In  classical  isolated  airfoil  and  axial  flow  cascade  unsteady  aerodynamic  models, 
upstream  airfoil  wakes  are  modeled  by  convected  harmonic  gusts.  Figure  5.4.  The 
behavior  of  these  convected  gusts  was  discussed  in  Chapter  4,  with  the  radial  and 
tangential  unsteady  velocity  components  shown  again  in  Equations  5.18  and  5.19. 

Vr  =  L^  =  ilf&L^2-l£  exp(A  In  r  +  Br2)  r1 '  k2  dr  - 


expfAlnr  +  Br2)/-1  +  k2dr  exp(ik(0 - r ))  + 


ji^L  rk2  - 1  +  rk2  - 1)  exp(ik(0 - 1 )) 


-  .  - .  roQo  rk2  - 1 

dr  2 


r 

I  exp(A  In 


r  +  Br2)r1*k2dr  + 


exp(A  In  r  +  B  r2)  r1  +  k2  dr  exp(ik(0- 1 ))  + 


('  rk2 ' 1  +  ^  r'k2  ‘ !)  - 1 » 
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From  these  velocity  components,  the  normal  velocity  induced  by  the  wakes  on  the  airfoil 
surfaces  is  computed.  Equation  5.20. 

v„  =  vrsin  Oo  -  VqCos  Oq  5.20 


The  expressions  for  the  radial  cascade  wake  behavior,  Equations  5.18  and  5.19,  cannot 
be  evaluated  for  k2  <  2,  due  to  the  fact  that  Qq  cannot  be  chosen  so  that  the  wake  velocities 

will  remain  bounded  as  r  -*».  Because  of  this,  the  standard  axial  cascade  upwash, 
Equation  5.24,  is  used  in  the  transformed  axial  cascade  at  low  reduced  frequency. 

v  =  v0  exp(ik'Ax  -  ikAt)  5 .24 


where  the  reduced  frequency,  kA,  is  equal  to  the  transformed  reduced  frequency  in 
Equation  5.14,  and  k’A  is  a  transformed  boundary  condition  reduced  frequency.  This 
separate  reduced  frequency  occurs  due  to  the  deceleration  of  the  steady  flow  with  radius, 
thus  decreasing  the  wavelength  of  the  wake. 

The  transformation  reduced  frequency.  Equation  5.14,  is  necessary  due  to  the  chord 
length  change  through  the  transformation  and  also  the  change  in  velocity  at  the  leading 
edge.  The  reduced  frequency  in  the  boundary  condition.  Equation  5.24,  is  the  ratio  of  the 
chord,  c,  to  the  gust  wavelength,  U/co,  thus  determining  the  number  of  wake  wavelengths 
on  the  vane.  Due  to  the  velocity  decrease  with  radius  in  the  radial  cascade,  the  gust 
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wavelength  becomes  smaller  toward  the  trailing  edge.  The  transformed  boundary  condition 
reduced  frequency  is  then 


—  Jc  (2  ~  At/T2) 

2(1  -  Ar/T2)  5.25 

The  boundary  condition,  Equation  5.20,  is  only  necessary  at  low  k,  with  the  more 
physically  realistic  wake  velocities  in  Equations  5.18  and  5.19  used  at  higher  values  of  k 
which  would  be  encountered  in  real  machines. 


5,4,  Model,  VqificaBOD 

To  validate  the  transformation  technique  and  radial  cascade  modeling,  the  radial  cascade 
unsteady  aerodynamic  lift  and  moment  are  predicted  for  the  special  case  of  Ar/r2— >0.  In 
this  case,  the  radial  cascade  flowfield  approaches  that  of  an  axial  flow  cascade.  The  radial 
cascade  parameters  are  Nv  -  300,  a  =  60°,  A rfa  =  0.01  and  rea  =  0,  as  defined  in  Figure 
5.2.  For  the  axial  cascade,  s/c  »  1.0,  o  =  60°,  and  x^  =  0.  The  number  of  airfoils  in  the 
radial  cascade  was  chosen  to  give  a  transformed  axial  cascade  solidity  of  1.0. 

Figures  5.5  through  5.10  show  the  unsteady  lift  due  to  translation,  the  unsteady 
moment  due  to  torsion  and  the  unsteady  lift  due  to  a  transverse  gust  for  both  the  baseline 
axial  and  radial  flow  cascades  over  a  range  of  reduced  frequencies  and  interblade  phase 
angles.  The  predicted  unsteady  aerodynamics  of  the  baseline  radial  cascade  approach  that 
of  the  axial  cascade,  thereby  verifying  the  transformation  and  the  radial  cascade  unsteady 
aerodynamic  modeling. 


5,5-Rg5uUs 

The  fundamental  differences  between  radial  and  axial  cascade  unsteady  aerodynamics 
are  considered  by  utilizing  this  model  to  predict  radial  and  corresponding  axial  cascade 
unsteady  aerodynamics.  In  particular,  the  unsteady  aerodynamic  translational  lift,  torsional 
moment,  and  convected  gust  lift  as  well  as  the  corresponding  unsteady  pressure 
distributions  are  analyzed.  The  baseline  radial  cascade  parameters  are  Nv  =  12,  a  =  60°, 
ATfo  =  0.25  and  rea  =  0.  The  baseline  geometry  differs  from  the  previous  cases  only  in 
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that  the  vane  radial  extent,  Ar/r 2,  is  greater,  thereby  demonstrating  the  effects  of  the  radial 
geometry. 

Figure  5.11  shows  the  unsteady  aerodynamic  translational  lift.  Cl,  for  the  baseline 
radial  cascade  as  a  function  of  the  reduced  frequency  and  interblade  phase  angle.  For  this 
translational  response,  the  unsteady  lift  will  amplify  the  motion  when  the  real  part  of  Cl  is 
positive,  indicating  a  possible  flutter  condition.  For  both  the  axial  and  radial  cascades,  all 
points  shown  have  negative  real  values  and  thus  are  stable.  Also,  the  stability  margins  for 
both  cascades  are  similar.  The  unsteady  lift  amplitude  for  the  radial  cascade  is  slightly 
higher  in  general  than  the  axial  flow  cascade. 

Figure  5.12  presents  the  unsteady  aerodynamic  torsional  moment.  Cm.  for  the  baseline 
radial  cascade  as  a  function  of  the  reduced  frequency  and  interblade  phase  angle.  In  this 
case,  the  imaginary  part  of  Cm  determines  stability,  with  a  stable  cascade  for  Cm  <  0.  The 
stability  is  similar  for  both  the  axial  and  radial  cascades,  with  the  zero  reduced  frequency 
results  unstable  for  both  cascades,  but  all  other  reduced  frequency  configurations  stable. 
As  for  the  translational  response,  the  interblade  phase  angle  has  a  lesser  effect  on  the 
unsteady  moment  for  the  radial  cascade  than  for  the  axial.  The  unsteady  moment 
amplitudes  are  greater  for  the  axial  cascade  at  low  reduced  frequency  and  become  roughly 
equal  at  higher  reduced  frequency. 

The  unsteady  aerodynamic  response  of  the  baseline  cascade  to  a  convected  transverse 
gust  as  a  function  c  f  Jie  reduced  frequency  and  interblade  phase  angle  is  shown  in  Figure 
5.13.  Since  the  reduced  frequency  is  low  (k  <  2),  the  standard  axial  cascade  upwash, 
Equation  4.2,  is  used  in  the  transformed  axial  cascade.  As  with  the  translational  and 
torsional  responses.  Figure  5.10  shows  that  the  interblade  phase  angle  has  a  lesser  effect  in 
the  radial  cascade  than  the  axial  cascade.  The  unsteady  lift  amplitude  is  smaller  for  the 
radial  cascade  than  the  axial 

Figures  5.14  and  5.15  show  the  unsteady  airfoil  pressure  difference  Cap  due  to 
translation  and  torsion  for  a  single  reduced  frequency/interblade  phase  angle  combination 
from  each  of  Figures  5.1 1  through  5.13.  The  unsteady  pressure  difference  is  presented  for 
the  radial  cascade  and  for  the  axial  cascade  conditions  of  Figures  5.5,  5.7  and  5.9.  The 
unsteady  pressure  is  greater  for  the  radial  cascade  undergoing  plunging  motions  which  is 
consistent  with  Figure  5.11.  The  unsteady  pressure  due  to  pitching  is  greater  at  the  airfoil 
leading  edge  for  the  radial  cascade  and  greater  at  the  trailing  edge  for  the  axial  cascade. 

Figure  5.16  shows  the  unsteady  airfoil  pressure  difference  due  to  a  convected  gust 
upwash  for  radial  and  axial  flow  cascades.  In  this  figure,  k  =  2  and  <$>(,  =  0°,  with  the 
transformed  cascade  upwash  given  by  Equation  5.20.  In  Figure  5.17,  k  =  4  and  the  radial 
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cascade  airfoil  pressure  difference  with  the  wake  model  upwash.  Equations  5.18  and  5.19, 
is  shown.  The  axial  cascade  unsteady  pressure  is  greater  than  the  radial  cascade,  consistent 
with  Figure  5.13.  In  Figure  5. 17,  the  axial  cascade  unsteady  pressure  is  greater  than  the 
facial  cascade  with  the  standard  upwash,  Equation  5.20,  but  is  less  than  the  radial  cascade 
with  the  wake  model  upwash,  Equations  5.18  and  5.19.  Using  these  more  physically 
realistic  wake  expressions  increases  the  predicted  radial  cascade  lift 

In  Chapter  6,  the  unsteady  diffuser  vane  pressure  distribution  measured  in  the  Purdue 
Centrifugal  Compressor  Facility  will  be  presented.  The  preceding  analysis  will  be 
compared  with  the  measured  data. 
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Figure  5.1  Radial  Cascade  Conformal  Transformation 
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Figure  5.2  Flat  Plate  Cascade  Transformation 
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Figure  5.4  Shed  Vortical  Gusts  in  Axial  and  Radial  Turbomachines 
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Figure  5.6  Baseline  Radial  Flow 
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Figure  5.9  Baseline  Axial  Flow  Cascade  Convected  Gust  Lift 
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Figure  5.12  Radial  Flow  Cascade  Torsional  Moment 
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CHAPTER 6 

EXPERIMENTAL  RESULTS 


A  set  of  experiments  was  performed  in  the  Purdue  Research  Centrifugal  Compressor  to 
investigate  the  overall  impeller  and  diffuser  performance,  impeller  and  diffuser  steady 
flowfield,  unsteady  radial  diffuser  wake  behavior,  unsteady  diffuser  vane  surface  static 
pressure  and  the  occurrence  of  rotating  stall  and  surge.  These  quantities  are  measured  for 
varying  flow  rate,  Q,  number  of  diffuser  vanes,  Nv,  diffuser  vane  stagger,  o,  diffuser 
vane  leading  edge/impeller  exit  radius  ratio,  rjAb*  and  diffuser  vane  row  detuning  using 
nonuniform  circumferential  spacing.  The  various  machine  geometric  conditions  are  shown 
in  die  following  table. 


Table  6.1  Research  Compressor  Geometric  Configurations 


riAo  *  1.10 
ri/r0  =  1.15 

ti/t0=  1.10 

riAb®  110 
ri/r0  =  1.10 
C409T.  detuned) 
nha  =  1.15 

rj/r0  =  1.10 
ri/r0=  1.15 

ri/rQ=  1.10 
ri/r0=  1.15 

ri/r0=  1.10 

At  each  of  die  diffuser  geometric  configurations  listed  in  Table  6. 1 ,  steady  data  are  acquired 
over  a  range  of  flow  coefficients,  typically  from  <D  =  0.16to  0.33.  Additionally,  unsteady 
airfoil  surface  pressure  data  are  acquired  at  0  =  0.17, 0.22  and  0.30.  Laser  velocimeter 
measurements  are  made  for  the  vaneless  configuration  atO= 0.17, 0.22  and  0.30. 
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6.L.CompBS5or  Performance 

The  compressor  performance  was  determined  for  each  of  the  geometric  conditions  in 
Table  6.1.  The  performance  is  presented  as  the  pressure  rise  coefficient,  'P,  ami 
efficiency,  T|,  versus  the  flow  coefficient,  <D.  For  the  vaneless  diffuser,  the  results  are 
presented  for  the  impeller  and  also  for  the  impeller/diffuser  stage.  For  the  vaned  diffuser, 
only  the  stage  results  are  presented.  The  high  flow  coefficient  on  the  performance  plots 
corresponds  to  the  completely  open  throttle  valve,  whereas  the  low  flow  coefficients 
correspond  to  stall  or  surge  (typically  O  »  0.16).  The  stall  and  surge  points  will  be 
discussed  more  fully  in  Section  6.6. 

Generally,  vaned  diffusers  are  utilized  to  improve  the  stage  performance.  The  vaned 
diffuser  itself  imposes  a  performance  loss  and  can  adversely  affect  the  overall  operating 
range  as  well  as  the  stall  and  surge  margins.  The  benefits  are  due  to  its  ability  to  favorably 
affect  the  aft  impeller  and  vaneless  space  flowfields,  thereby  improving  overall  stage 
performance  and  sometimes  extending  the  machine  operating  range  into  either  lower  or 
higher  flow  rates. 

The  vaneless  diffuser  machine  performance  is  presented  in  Figure  6. 1 .  The  maximum 
pressure  rise  occurs  at  <X>  =  0.19,  with  the  maximum  efficiency  at  <D  =  0.22.  These  values 
are  also  typical  of  the  vaned  diffuser  performance.  There  is  typically  a  5%  loss  in  both  the 
total  pressure  and  efficiency  through  the  vaneless  diffuser  at  each  operating  point 

Figures  6.2  and  6.3  present  the  stage  performance  for  various  values  of  the  diffuser 
vane  stagger  angle,  O.  In  these  figures,  rj/r0  =  1.10,  and  Nv  *  30  and  15  respectively. 
For  Nv  =  30,  the  stage  performance  with  the  vaneless  diffuser  is  typically  similar  or  better 
than  that  with  the  vaned  diffuser.  The  stage  performance  at  a  =  70°  is  particularly  worse 
due  to  the  Urge  blockage  (vane  row  open  area/total  area).  For  the  vaned  diffuser  with 
Nv  =  15,  the  stage  performance  is  greater  than  for  the  vaneless  stage.  For  all  values  of  a, 
the  vaned  diffuser  yields  a  stage  performance  improvement  over  that  with  the  vaneless 
diffuser  for  <b  >  0.27.  At  the  maximum  efficiency  point,  there  is  a  3%  efficiency  gain  fc  r 
a  =  70°  over  the  vaneless  case.  The  performance  is  greater  at  higher  d>  for  smaller  values 
of  a  and  greater  at  lower  G>  for  larger  o.  This  is  due  in  part  to  the  incidence  angle  variation 
on  the  diffuser  vane  row  with  flow  rate.  Since  the  impeller  is  backswept,  the  absolute  flow 
angle,  Oq,  decreases  with  increasing  flow  rate,  increasing  the  incidence  on  the  diffuser 
vanes.  Therefore,  smaller  values  of  o  yield  improved  performance  at  higher  flow  rates 
with  higher  a  better  at  low  flow,  In  Figure  6.4,  ri/rQ  =  1.15  and  Nv  =  30.  The 
performance  is  similar  to  n/r0  =1.10  except  for  a  =  60°  where  the  performance  is 
significantly  improved.  This  behavior  is  not  understood  since  the  performance  at  a  =  50° 
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and  70°  is  slightly  greater  at  ri/r0  =  1.10  and  typically,  the  performance  is  improved  for 
closer  vane/impeller  spacing. 

Figures  6.5  through  6.7  present  the  stage  performance  for  various  numbers  of  diffuser 
vanes,  Nv.  In  Figure  6.5,  the  performance  is  slightly  greater  for  Nv  =  15,  except  for 
4»0.30,  where  presumably  the  greater  number  of  vanes  delays  separation  at  negative 
incidence  angles.  In  Figure  6.6,  the  performance  is  similar  for  Nv  =  30  and  15  but 
improved  at  higher  flow  rates  for  Nv  =  10.  In  Figure  6.7,  the  performance  improves  for 
Nv  =  15. 

As  the  vane  stagger  angle,  a,  is  increased,  the  area  blocked  by  the  diffuser  vane  row 
increases  with  the  flow  area  between  die  vanes  decreasing.  This  reduces  the  flow  area  with 
respect  to  die  surface  area,  thus  increasing  die  relative  effects  of  skin  friction.  Therefore,  at 
lower  values  of  a,  the  number  of  diffuser  vanes  has  a  lesser  effect  but  for  higher  o,  the 
blocked  area  increases  with  more  vanes  yielding  a  performance  penalty.  If  the  number  of 
vanes  is  small,  the  skin  friction  is  less  but  the  ability  of  the  vane  row  to  favorably  affect  the 
flow  and  to  avoid  separation  is  diminished.  Therefore,  competing  effects  due  to  the 
number  of  vanes  and  vane  stagger  angle  exist  which  must  be  determined  from  the  desired 
performance  and  operating  range. 

In  Figure  6.8,  the  effect  of  radius  ratio,  rifto,  is  shown  for  o  =  60°.  A  significant  gain 
in  total  pressure  and  efficiency  is  seen  for  ri/r0  =  1.15  as  compared  to  ri/r0  =  1.10  or  the 
vaneless  case,  with  a  maximum  efficiency  gain  of  3%.  For  o  =  50°  and  a  =  70°,  the 
change  in  performance  due  to  radius  ratio  is  minimal. 

The  performance  for  the  detuned  diffuser  vane  row  using  nonuniform  circumferential 
spacing  is  also  determined  To  accomplish  this,  every  other  blade  mounting  cam  is  rotated 
in  the  opposite  direction  of  its  neighbor,  yielding  a  nonuniform  circumferential  airfoil 
spacing.  Figure  6.9.  An  even  number  of  blades  must  be  used  and  only  one  detuning  ratio 
is  possible  for  each  o  and  riAo.  tine  to  the  cam  geometry.  Figures  6.10  and  6.1 1  show  the 
effect  of  detuning  far  a  =  60°  and  rj/r0  =  1.10.  For  Figure  6.10,  Nv  =  30  with  40% 
detuning,  indicating  every  other  airfoil  has  been  angularly  displaced  40%  of  the  uniformly 
spaced  cascade  angle.  For  this  case,  the  performance  change  is  minimal,  with  all 
differences  within  the  experimental  uncertainty.  For  Figure  6.11,  Nv  10  with  13% 
detuning.  A  slight  increase  in  the  efficiency  and  pressure  rise  coefficient  is  seen  for  lower 
flow  rates.  For  the  tuned  vane  row  with  o  =*  60°  and  n/r0  *  1. 10,  the  performance  with 
Nv  -  10  is  decreased  compared  to  that  for  higher  numbers  of  vanes  at  lower  flow  rates. 
Perhaps  the  detuning  yields  an  effect  similar  to  increasing  Nv.  since  every  other  blade  is 
more  closely  spaced. 
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In  summary,  the  vaned  diffuser  can  improve  the  stage  performance  over  that  of  a 
vaneless  diffuser.  The  performance  improvement  is  greatest  when  the  diffuser  vane 
incidence  is  low  and  the  blockage  is  small.  The  diffuser  vane  stagger  angle  and  the  number 
of  vanes  were  the  most  important  parameters. 


6.2  Impeller  Steady  Loading 

Impeller  blade  static  pressure  measurements  are  made  to  determine  the  steady  loading 
on  the  rotating  impeller.  This  is  accomplished  using  the  rotating  Scanivalve  described  in 
Chapter  3.  Measurement  locations  are  shown  in  Figure  2.5.  The  pressure  and  suction 
surface  static  pressure  coefficient,  Cp,  and  the  pressure  difference  coefficient,  CAP  are 
plotted  versus  the  nondimensional  impeller  radius,  r/rQ  along  three  streamlines. 

Figures  6.12  through  6.14  show  the  impeller  blade  surface  static  pressure  distributions 
for  ®  =  0.17,  0.22  and  0.30  for  the  vaneless  diffuser.  At  ®  =  0.17,  the  blade  loading 
indicated  by  the  pressure  difference  coefficient,  CAP,  is  uniform  at  the  hub  and  mean  with 
the  tip  having  higher  loading  at  impeller  inlet  and  exit.  The  higher  loading  at  the  impeller 
inlet  tip  is  likely  due  to  positive  blade  incidence  caused  by  the  low  axial  velocity.  The  blade 
loading  at  the  maximum  efficiency  point,  ®  =  0.22,  is  uniform  for  each  streamline  with  the 
the  loading  slightly  less  than  that  at  ®  =  0.17  along  the  mean  and  tip  streamline.  The 
impeller  blade  loading  at  ®  =  0.30  is  less  uniform  than  the  previous  cases,  although  higher 
in  magnitude.  At  the  impeller  inlet,  a  low  pressure  on  the  pressure  surface  and  higher 
pressure  on  the  suction  surface  occurs,  especially  toward  the  blade  tip,  indicating  a  slight 
negative  incidence  separation.  Generally,  the  blade  loading  increases  with  flow  rate,  due  to 
the  increased  coriolis  acceleration. 

In  Figure  6.15,  the  impeller  blade  surface  static  pressure  and  pressure  difference  for 
r|/r0  *  1.10,  a  =  50°,  Nv  =  30  and  4>  =  0.30  are  presented  along  with  the  vaneless  data  at 
®  =  0.30.  The  impeller  blade  loading  indicated  by  is  greater  for  the  vaneless  diffuser 
over  the  entire  blade  except  at  the  impeller  inlet  tip,  where  the  loading  with  the  vaned 
diffuser  is  slightly  greater.  Although  the  loading  is  higher  with  the  vaneless  diffuser,  the 
overall  static  pressure,  Cp,  is  generally  greater  for  the  vaned  diffuser.  In  Figure  6. 16,  the 
impeller  blade  surface  static  pressure  and  loading  for  rj/r0  =  1.10,  a  =  60°,  Nv  =  30  and 
$  =  0.22  is  presented  along  with  the  vaneless  data.  The  vaned  diffuser  has  caused  a  slight 
increase  in  impeller  loading  at  the  tip  and  mean  blade  positions.  The  results  for  c  =  70°  are 
similar. 
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For  Figures  6. IS  and  6.16,  the  diffuser  vane  incidence  angle  is  low  and  therefore  the 
diffuser  vane  loading  indicated  by  the  steady  pressure  difference  is  small.  For  O  =  0.30  in 
Figure  6.15,  the  performance  indicated  by  the  efficiency  and  pressure  rise  coefficients  is 
higher  for  the  vaned  diffuser.  For  Figure  6. 15  with  d>  =  0.22,  the  performance  is  the  same 
for  the  vaned  and  vaneless  diffusers.  At  all  other  flow  rates  and  diffuser  geometries,  no 
effect  of  diffuser  geometry  is  seen  on  the  impeller  blade  static  pressure  or  loading. 

Therefore,  the  impeller  loading  indicated  by  the  blade  surface  static  pressure  differential 
generally  increases  with  increased  flow  rate.  The  impeller  blade  loading  is  only  affected 
when  30  diffuser  vanes  are  used  and  the  diffuser  vane  incidence  is  low  producing  low  vane 
loading. 


6.3  Diffuser  Vane  Steady  Loading 

Diffuser  vane  surface  static  pressure  measurements  are  made  at  the  vane  mean  axial 
position  and  chordwise  locations  indicated  in  Figure  3.2.  The  pressure  and  suction  surface 
steady  pressure  coefficients,  Cp,  and  the  pressure  difference  coefficient,  CAP  are  presented 
as  a  function  of  nondimensional  chordwise  distance. 

In  Figure  6.17,  the  diffuser  vane  static  pressure  is  presented  for  x\h0  =  1.10,  c  =  60° 
and  Nv  =  15.  Three  curves  are  shown  on  each  plot,  each  for  a  separate  value  of  O  and 
therefore  different  diffuser  vane  mean  incidence  angles,  i.  At  d>  =  0.30,  the  incidence  is 
slightly  positive,  i  =  4.5°,  and  the  vane  loading  is  small,  Cl  =  -0.146.  As  is  expected, 
decreasing  the  incidence  to  greater  negative  angles  cause  the  loading  to  increase  negatively. 
The  vane  loading  for  a  given  Nv  is  independent  of  c  and  therefore  not  shown  in  the 
figures. 

In  Figure  6.18,  ri/r0  =  1.10,  a  =  60°,  O  =  0.22  and  i  =  -5.8°.  Three  curves  are 
presented  per  plot  for  Nv  =  30, 15  and  10  vanes.  Decreasing  the  vane  number  causes  the 
loading  to  be  more  negative.  This  behavior  is  typical  for  other  diffuser  geometries 
investigated. 

The  effect  of  detuning  on  the  diffuser  vane  row  is  shown  in  Figure  6.19,  with  a  40% 
detuned  vane  row  as  defined  in  Figure  6.9.  In  this  figure,  ri/r0  =  1.10,  a  =  60°,  Nv  =  30, 
d>  =  0.30  and  i  =  4.5°.  The  loading  on  the  detuned  vanes  is  changed  slightly,  with  the 
leading  edge  pressure  surface  pressure  affected  most  significantly.  The  loading  on  vane  1 
has  increased  above  the  tuned  vaned  whi<e  the  loading  on  vane  2  has  decreased.  At  lower 
flow  rates,  the  pressure  change  due  to  detuning  is  less  significant.  The  case  with  13% 
detuning  and  Nv  =  10  is  not  shown,  since  no  change  in  loading  due  to  detuning  occurs. 
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It  is  interesting  to  note  that  in  cases  where  the  greatest  performance  increase  is  realized 
utilizing  the  vaned  diffuser,  the  vanes  are  typically  lighdy  loaded.  This  lends  credibility  to 
the  belief  that  the  vaned  diffuser  increases  performance  by  affecting  the  aft  impeller  and 
vaneless  space  flewfields  (Fisher  and  Inoue,  1981),  thereby  reducing  secondary  flows  and 
separation.  Thus,  if  a  vaned  diffuser  is  utilized  to  increase  efficiency,  the  vane  stagger 
should  be  equal  to  the  absolute  flow  angle  at  the  operating  point,  minimizing  the  loading. 

In  summary  ,  the  diffuser  vane  loading  indicated  by  the  surface  static  pressure  and 
steady  lift  coefficient,  is  mainly  dependent  on  the  vane  incidence  and  number  of  vanes. 
Changing  the  stagger  angle  has  a  minimal  effect  on  the  loading  if  the  incidence  is  held 
constant 

6.4  Hadial  Diffuser  Velocity  Field 

LDV  velocity  measurements  are  made  at  each  of  the  axial  and  radial  positions  shown  in 
Figure  6.20  for  the  vaneless  diffuser  at  0  =  0.17,  0.22  and  0.30.  The  instantaneous 
velocity  is  measured  and  the  mean  radial  and  tangential  velocity  in  each  measurement  bin 
around  the  impeller  circumference  is  determined,  as  described  in  Section  3.3  and  depicted 
tn  Figures  6.21  and  6.22. 

Figures  6.23  through  6.25  show  the  steady  time-averaged  radial  and  tangential 
velocities,  Ur  and  Ue,  nondimensionalized  by  the  diffuser  inlet  velocities,  Uri  and  Uej. 
These  are  presented  at  each  axial  position  as  a  function  of  nondimensional  radius,  r  =  r/ri', 
where  rj'  is  the  vaneless  diffuser  inlet  depicted  in  Figure  6.22.  Each  figure  represents  a 
separate  flow  coefficient  and  the  incompressible,  inviscid,  ideal  solution,  Ur/Uri  =  1/r  and 
Ue/Uei  =  1/r  are  also  shown. 

For  0  =  0.22  and  0.30,  the  tangential  velocities  follow  the  ideal  solution,  whereas  the 
radial  velocities  tend  to  remain  constant,  except  at  tne  tip  for  0  =  0.30,  where  a  decrease  is 
seen.  This  behavior  is  presumably  due  to  boundary  layer  formation  on  the  diffuser  walls, 
decreasing  the  velocity  near  the  wall  and  increasing  the  core  velocity  to  satisfy  continuity. 
Viscous  effects  and  mixing  appear  to  have  little  effect  on  the  tangential  velocity  since  it 
tends  to  follow  the  ideal  behavior.  Therefore,  viscous  effects  in  the  diffuser  appear  to  be 
confined  to  the  wall  region  and  have  not  penetrated  to  the  measurement  locations. 

For  0  =  0.17  in  Figure  6.25,  the  flow  shows  increased  departure  from  the  ideal 
solution.  This  operating  point  is  at  incipient  stall  which  is  discussed  in  Section  6.6 
(0  =  0.16  for  the  vaneless  diffuser).  Therefore,  turbulent  mixing  and  separation  are  more 
likely.  A  slight  decrease  in  the  tangential  velocity  from  the  ideal  solution  is  observed.  The 
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radial  velocity  shows  a  significant  departure  with  an  increase  at  the  tip  and  mean  positions 
and  a  decrease  at  the  hub.  At  the  two  largest  measurement  radii,  a  more  pronounced 
departure  from  the  theory  is  seen  in  the  radial  and  tangential  velocities  indicating  a  possible 
diffuser  wall  separation. 

In  Figures  6.26  through  6.34,  a  vector  plot  of  the  vaneless  diffuser  total  velocity, 
ensemble  averaged  using  the  data  from  the  individual  impeller  passages,  is  shown  for  each 
flow  coefficient  and  axial  position.  This  geometry  of  these  figures  is  depicted  in  Figure 
6.22.  The  wheel  speed  magnitude  and  direction  is  shown  for  comparison  and  the  scale  is 
the  same  for  each  figure.  In  Figures  6.26  through  6.28,  the  results  for  O  =  0.30  are 
presented.  The  circumferential  nonuniformity  of  the  flow  can  be  observed,  especially  at 
r  »  1.0  and  1.06  at  the  tip  and  mean  positions.  The  impeller  exit  flow  angle,  oto,  is 
greatest  at  the  tip  due  to  the  greater  impeller  tip  radius  and  thus  greater  wheel  speed.  This 
behavior  is  consistent  at  each  radius.  The  velocity  magnitude  shows  the  opposite  behavior 
with  the  hub  velocities  being  the  greatest  Similar  behavior  is  observed  in  Figures  6.29 
through  6.31  for  $  =  0.22,  although  the  velocity  magnitude  is  more  uniform  from  tip  to 
hub.  In  Figures  6.32  through  6.34  for  $  =  0.17,  the  circumferentially  nonuniform  flow  is 
less  pronounced  than  in  the  previous  figures.  The  flow  angle  is  greatest  at  the  tip  for 
r  =  1.0  and  1.06,  but  this  behavior  is  reversed  at  larger  radii.  At  the  hub  position,  a  sharp 
velocity  decrease  is  seen  at  r  =  1.32  and  1.45,  which  is  also  observed  in  Figure  6.25.  It 
appears  that  a  thick  boundary  layer  has  developed  or  that  separation  is  occurring  on  the 
diffuser  hub  wall. 

In  Figures  6.35  through  6.43,  the  vaneless  diffuser  unsteady,  circumferentially 
nonuniform  velocity  is  presented  by  subtracting  the  steady  time  averaged  velocity  from 
each  vector  in  the  previous  figures.  These  figures  more  clearly  show  the  effect  of  flow  rate 
and  radius  on  the  unsteady  radial  diffuser  wake  behavior.  In  each  figure,  the  steady  radial 
velocity  magnitude  and  direction  at  the  same  axial  location  is  shown  for  comparison.  At  the 
mean  axial  positior  and  r  =  1.0  for  each  d>,  Figures  6.36,  6.39  and  6.42,  the  first 
harmonic  circumferential  impeller  exit/diffuser  entrance  unsteady  static  pressure  is  shown 
scaled  by  pUi.  Since  the  unsteady  vaned  diffuser  airfoil  surface  static  pressures  are 
nondimtnsionalized  by  pUiv„,  the  relative  effect  of  the  unsteady  impeller  exit/diffuser 
entrance  static  pressure  on  the  vaned  diffuser  can  be  estimated.  The  plot  scales  for 
4>  =  0.22  and  =  0.17  are  2x  and  4x  respectively  that  of  O  =  0.30. 

In  Figures  6.35  through  6.37  for  O  =  0.30,  the  wake  is  mainly  confined  to  the  tip  and 
mean  positions.  It  can  be  seen  that  the  wake  has  nearly  disappeared  at  r  =  1.45,  which  is 
typical  for  all  of  the  measurements  made.  For  O  =  0.22  in  figures  6.38  through  6.40,  the 
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wake  is  also  confined  mainly  to  the  tip  and  mean.  The  wake  velocity  vectors  are  less 
parallel  at  r  =1.0  than  for  =  0.30.  At  0  =  0.17  in  Figures  6.41  through  6.43,  the  wake 
velocity  vectors  tend  to  be  more  aligned  with  the  absolute  flow  whereas  the  vectors  at 
O  =  0.30  and  0.22  are  more  aligned  with  die  flow  relative  to  the  impeller. 

A  comparison  of  the  vaneless  radial  diffuser  unsteady  velocities  shown  in  Figures  6.35 
through  6.43  is  made  with  the  theory  developed  in  Chapter  4.  This  is  accomplished  by 
taking  the  fir  ♦  harmonic  unsteady  radial  and  tangential  velocities,  scaled  by  the  magnitude 
at  the  van  dess  diffuser  inlet,  |v  |  ={  | ViP  +  fveP)1/  and  plotting  the  magnitude  and  phase 
along  with  theory.  Because  measurements  were  not  made  immediately  at  the  impeller  exit, 
the  vaneless  radial  diffuser  inlet,  ri\  will  be  considered  the  impeller  exit  Before  these 
results  are  presented,  it  is  useful  to  consider  the  parameters  which  describe  wakes  in 
classical  linear  airfoil  theory  for  axial  flow  cascades. 

Typically,  wakes  shed  from  turbomachine  blade  rows  are  considered  small  amplitude 
sinusoidal  gusts  superimposed  on  the  steady  relative  flow  leaving  the  rotor,  Figures  4.1 
and  6.44.  The  upper  example  in  Figure  6.44  depicts  the  classical  linear  theory  gust  where 
the  velocity  vectors  are  parallel  and  aligned  with  the  relative  flow  angle,  Po-  In  this  type  of 
gust,  the  gust  phase  angle  or  phase  angle  between  vr  and  v0  is  <j>w  =  180°,  and  the  ratio, 
vJvq  =  tan  po.  In  the  bottom  left  example,  <J>W  =  180°  but  v,/v0  #  tan  p0.  Here,  the 
velocity  vectors  are  parallel  but  not  aligned  with  the  relative  flow  angle,  p0.  In  the  bottom 
right  example,  Vv0  =  tan  p0,  but  <j>w*180°.  In  this  example,  the  velocity  vectors  are  no 
longer  parallel. 

Figures  6.45  through  6.47  show  the  mean  axial  position  first  harmonic  unsteady  radial 
and  tangential  velocity  components  compared  with  the  theory  from  Chapter  4.  In  the 
theory  as  mentioned  previously,  <|>w  =  180°  and  vr/v0  =  tan  p0,  although  the  theory  can 
accommodate  <j>w  *  180°  and  v,/v0  *  tan  p0.  For  the  experimental  data,  the  velocity  ratio 
and  gust  phase  angle  as  well  as  the  absolute  and  relative  flow  angles,  cto  and  p0,  are 
calculated  at  the  diffuser  inlet,  rj’  =  40.4  cm  (15.9  in.).  In  Figure  6.45,  <I>  =  0.30,  the  gust 
phase  angle,  <J>W,  is  160°  and  v,/v0  =  tan  p0.  In  a  general  sense,  the  magnitude  of  vr 
follows  the  theoiy.  However,  the  magnitude  of  vq  decreases  more  rapidly  with  r  than  the 
theory.  Thephaseof  both  vr  and  ve  vary  less  rapidly  with  r  than  in  the  theory.  In  Figure 
6.46,  O  =  0.22,  the  gust  phase  angle,  <|>w,  is  148°  and  again,  v,/v0  =  tan  P0.  As  in  the 
previous  figure,  the  radial  velocity  magnitude  generally  follows  the  theory  with  the 
tangential  velocity  decreasing  more  rapidly  with  r.  The  phase  varies  less  rapidly  with  r  for 
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the  data  than  for  the  theory.  In  Figure  6.47,  <I>  =  0.17,  the  gust  phase  angle,  <|>w»  is  60° 
and  Vve  » tan  f}0.  The  radial  and  tangential  velocity  magnitudes  decrease  more  rapidly 
v/ith  r  than  the  theory  with  poor  phase  correlation. 

It  can  be  seen  in  Figures  6.45  through  6.47  that  the  correlation  of  the  data  with  the 
theory  deteriorates  as  d>  decreases  and  the  gust  phase  angle,  $w,  departs  from  180°. 
Therefore,  as  the  wake  velocity  departs  from  the  linear  theory  profile  in  Figure  6.44,  the 
correlation  with  the  theory  becomes  worse. 

In  Figure  6.48,  Po  =  -53°  and  the  theoretical  <j>w  is  -160°,  equal  to  the  data  phase  angle. 
This  is  accomplished  by  introducing  a  phase  angle  between  Vn>  and  veo  in  the  boundary 
condition.  Equations  4.5.  The  phase  correlation  is  somewhat  improved  due  to  the  fact  that 
the  impeller  exit  relative  flow  angle  is  probably  closer  to  the  backsweep  angle  of  51°. 
However  the  amplitude  correlation  is  unimproved,  with  an  oscillation  in  the  theory 
amplitude  at  lower  r. 

It  is  hypothesized  that  the  presence  of  the  circumferential  unsteady  static  pressure 
gradient  plays  an  important  role  in  wake  behavior.  Because  an  unsteady  circumferential 
static  pressure  gradient  is  generated  by  potential  effects  due  to  impeller  rotation.  Section 
1.2,  with  the  velocity  nonuniformity  generated  by  viscous  effects,  the  wake  analysis  does 
not  directly  account  for  it.  Henderson  (1991)  has  shown  unsteady  circumferential  pressure 
gradients  to  exist  in  axial  machines  significant  distances  downstream  of  the  rotor,  and  that 
wake  behavior  and  downstream  airfoil  forced  response  are  significantly  different  than  in 
classical  linear  theory.  In  their  experiment,  linear  theory  type  gusts  are  only  generated  with 
perforated  plates  which  have  minimal  lift  and  no  large  areas  of  separation.  Therefore, 
surfaces  producing  lift  such  as  exist  in  any  turbomachine  will  produce  circumferential 
pressure  gradients  that  affect  wake  behavior  and  downstream  airfoil  response. 

In  summary,  it  is  seen  that  the  steady  circumferential  velocity  generally  follows  the 
ideal  theory  whereas  an  acceleration  of  the  steady  radial  velocity  due  to  boundary  layer 
thickening  causes  departure  from  theory.  The  unsteady  wake  velocities  decrease  rapidly 
with  radius,  with  the  wake  velocities  disappearing  at  a  nondimensional  radius  of  1.45.  The 
measured  unsteady  wake  velocities  generally  decrease  more  rapidly  with  increasing  radius 
than  predicted  by  the  theory.  The  phase  angle  between  the  unsteady  radial  and  tangential 
velocities  and  the  unsteady  circumferential  pressure  gradient  are  hypothesized  to  be 
important  factors  in  the  departure  from  theory. 
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6.5  Diffuser  Vane  Unsteady  Loading 

Unsteady  diffuser  vane  static  pressure  measurements  are  made  at  each  of  the  geometric 
configurations  listed  in  Table  6.1  at  <D  =  0.17,  0.22  and  0.30.  The  unsteady  lift 
coefficients  Cl  and  Cls»  defined  in  Equations  3.14  and  3.15,  are  presented  along  with  the 
steady  lift  coefficient.  Cl,  to  evaluate  the  diffuser  vane  loading  as  a  function  of  diffuser 
geometry  and  flow  rate.  The  diffuser  vane  surface  static  pressure  coefficients  and  pressure 
difference  coefficients  are  presented  along  with  the  theory  developed  in  Chapter  5. 

Figures  6.49  through  6.54  present  the  unsteady  and  steady  lift  coefficients,  Cls  and 
Cl,  as  a  function  of  flow  rate  for  the  various  diffuser  geometries.  The  unsteady  lift 
coefficient  Cls  is  the  Erst  harmonic  unsteady  lift  amplitude  nondimensionalized  by  the 
dynamic  pressure  based  on  wheel  speed,  Equation  3.15.  This  yields  an  absolute  measure 
of  the  vane  loading  as  a  function  of  the  machine  operating  point 

In  Figure  6.49,  the  unsteady  and  steady  lift  for  ri/r0  =  1.10  and  Nv  =  30  is  presented 
for  the  various  stagger  angles.  The  unsteady  lift  increases  markedly  with  G>  and  with 
increasing  stagger,  a  The  steady  lift  changes  more  rapidly  with  $  as  a  is  increased,  with 
Cl  becoming  more  positive  with  increasing  d>  and  a 

In  Figure  6.50,  the  radius  ratio  has  been  increased  to  ri/r0  =  1.15.  The  unsteady  lift 
changes  less  rapidly  with  d>  and  a,  although  the  lift  atd>  =  0.30  is  greater  for  cr  =  60°  and 
70°  in  Figure  6.50  than  for  Figure  6.49.  For  a = 60°,  the  unsteady  lift  is  greater  in  Figure 
6.50  with  ri/rQ  =  1.15  for  all  <D.  The  increased  lift  at  higher  ri/r0  is  surprising  since  the 
wake  velocities  are  decreasing  with  increased  radius.  The  steady  lift  behavior  is  similar  in 
both  cases.  It  is  interesting  to  note  that  the  steady  lift  for  a  =  50°  attains  a  minimum  at 
4>  =  0.22  and  further  decrease  in  <J>  and  thus  the  negative  incidence  has  minimal  effect  on 
the  lift  It  is  likely  that  the  vaned  diffuser  begins  to  stall  at  this  point.  This  behavior  is 
observed  for  all  diffuser  geometries. 

Figure  6.51  is  similar  to  Figure  6.49  but  with  Nv  =  15.  The  reduction  in  vane  number 
has  a  significant  effect  on  the  unsteady  lift  behavior.  The  unsteady  lift  is  typically  greater 
than  for  Nv  =  30,  however  the  effect  of  0  on  the  lift  is  different  at  each  value  of  <J.  For 
a  =  50°,  the  unsteady  lift  is  greatest  at  high  d>,  for  a  =  60°,  the  unsteady  lift  is  greatest  at 
moderate  O.and  for  c  =  70°,  the  unsteady  lift  is  greatest  at  low  0.  Each  of  these  stagger 
and  flow  rate  combinations  corresponds  closely  to  the  point  of  minimum  steady  incidence. 
The  steady  lift  is  shifted  in  the  negative  direction  but  shows  similar  behavior  to  the  previous 
figures.  In  Figure  6.52,  the  unsteady  and  steady  lift  is  shown  for  ri/rQ  =  1.15  and 
Nv  =  15.  Only  measurements  for  a  =  60°  are  made.  The  unsteady  lift  is  significantly 
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reduced  by  increasing  ri/r0  to  1.15.  At  this  radius  ratio,  the  unsteady  lift  is  a  minimum  for 
minimum  incidence,  as  opposed  to  a  maximum  as  in  the  previous  figure. 

In  Figures  6.53  and  6.54,  the  effect  of  detuning  on  the  unsteady  lift  is  presented  for 
Nv  =  30  with  40%  detuning,  and  Nv  =  10  with  13%  detuning  respectively.  For  all 
measurement  points,  the  unsteady  lift  is  increased  on  one  or  both  of  the  vanes.  It  cannot  be 
concluded  from  these  cases  that  detuning  will  always  produce  an  adverse  effect  on  the  lift, 
although  in  general  it  is  not  always  beneficial.  The  detuning  affects  the  steady  lift  less 
significantly  than  does  the  incidence. 

A  comparison  of  the  unsteady  airfoil  surface  pressure  data  is  made  with  the  unsteady 
radial  cascade  airfoil  theory  developed  in  Chapter  5.  The  experimental  diffuser  vane 
unsteady  static  pressure  coefficients,  Cp,  are  presented  for  both  the  suction  and  pressure 
surfaces,  along  with  the  experimental  and  theoretical  pressure  difference  coefficient,  Cap, 
as  defined  in  Chapter  3.  Pressure  coefficients  for  ri/rD  =  1.15  have  been 
nondimensionalized  by  the  transverse  velocity  at  n/r0  =  1.10  in  order  to  evaluate  the  radius 
ratio  effect  and  since  no  velocity  measurements  are  made  at  ri/r0  =  1.15.  Since  the  phase  is 
different  at  the  blade  leading  edge  than  at  the  velocity  measurement  location,  a  phase 
correction  is  made.  The  time  taken  for  the  wave  to  be  convected  from  the  measurement 
point  to  the  vane  leading  edge  is  computed  using  the  steady  velocity.corrected  for  phase. 

Before  malting  a  comparison  with  the  theory,  it  is  useful  to  consider  the  wake  velocity 
profiles  given  in  the  vector  plots  of  Figures  6.36, 6.39  and  6.42  along  with  the  data/theory 
comparison  in  Figures  6.45  through  6.47.  In  Figures  6.36  and  6.39  for  0  =  0.30  and 
0.22,  the  wake  velocity  vectors  are  nearly  parallel  with  the  steady  relative  flow,  and 
therefore  nearly  transverse  to  the  airfoils  in  the  radial  diffuser.  In  Figure  6.42  for 
0  =  0.17,  the  wake  velocity  vectors  are  nearly  parallel  with  the  absolute  steady  flow,  and 
therefore  nearly  parallel  to  the  airfoil  chord,  with  the  transverse  unsteady  velocity 
component  very  small.  It  is  observed  in  the  data/theory  comparison  of  Figures  6.45 
through  6.47  that  the  tangential  velocity  component  decreases  with  radius  more  rapidly  than 
predicted  by  the  theory.  However,  the  radial  component  showed  reasonable  magnitude 
correlation  in  a  general  sense  for  0  =  0.30  and  0.22.  Since  the  radial  velocity  is  the  major 
component  in  the  airfoil  transverse  velocity,  especially  at  higher  a,  the  deficiencies  in  the 
radial  diffuser  wake  theory  will  not  necessarily  be  a  major  source  of  error  in  the  unsteady 
airfoil  theory. 

Figures  6.55  through  6.65  present  the  unsteady  first  harmonic  airfoil  surface  static 
pressure  distributions.  On  the  lefthand  axes  of  each  figure,  the  amplitude  and  phase  of  the 
unsteady  static  pressure  coefficient,  Cp,  is  plotted  versus  nondimensional  vane  chord  for 
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pressure  and  suction  surfaces.  On  the  righthand  axes,  the  amplitude  and  phase  of  the 
unsteady  pressure  difference  coefficient.  Cap,  is  presented  along  with  the  theory  versus 
nondimensional  chord.  The  legend  presents  the  experimental  and  theoretical  lift 
coefficients,  Cl  and  CLt»  as  defined  in  Equation  3.14,  along  with  the  lift  coefficient,  Cls, 
defined  in  Equation  3.15. 

Figures  6.55  through  6.57  show  the  first  harmonic  unsteady  static  pressure  distribution 
for  O  =  0.30,  0.22  and  0.17,  respectively,  with  Nv  =  30,  ri/r0  =  1.10  and  a  =  60°.  For 
d>  =  0.30,  the  amplitude  of  Cp  is  high  at  the  leading  edge  for  pressure  and  suction  surfaces 
as  is  typical  in  steady  and  unsteady  cascade  flows.  The  pressure  surface  phase  is  nearly 
constant  along  the  airfoil  c>,  rd  with  some  variation  in  the  suction  surface  phase.  Good 
correlation  is  seen  in  the  phase  of  Cap  with  the  theory.  However,  the  theory  generally 
oveipredicts  the  amplitude  of  Cap  by  a  factor  of  2  or  more. 

In  Figure  6.56  for  <B  =  0.22,  the  amplitude  of  Cp  is  somewhat  higher  than  in  the 
previous  figure  but  the  behavior  is  similar.  The  phase  behavior  of  Cp  is  similar  for  the 
suction  surface  but  the  pressure  surface  shows  more  variation.  The  phase  correlation  of 
Cap  has  deteriorated  with  a  greater  phase  decrease  in  the  data  toward  the  trailing  edge.  The 
steady  flow  incidence,  i,  is  negative  for  this  case  and  positive  for  the  previous  case. 
Greater  phase  decrease  at  negative  incidence  is  seen  for  all  measurements  made  with 
Nv  =  30.  The  amplitude  is  closer  to  the  theory  at  the  airfoil  trailing  edge  than  in  the 
previous  case.  However,  the  theory  still  generally  overpredicts  the  amplitude  by  a  factor  of 
2  or  more. 

Figure  6.57  presents  the  unsteady  pressure  for  =  0.17.  The  amplitude  and  phase 
behavior  of  Cp  is  similar  to  the  previous  figure  but  the  amplitude  is  much  greater.  The 
phase  decrease  in  the  data  for  Cap  occurs  due  to  negative  incidence.  The  amplitude  of  Cap 
is  now  underpredicted  by  the  theory.  This  behavior  is  generally  seen  for  all  geometries  at 
<I>  =  0.17,  with  the  data  greater  than  the  theory  over  most  of  the  airfoil  chord.  This  is 
attributed  to  the  small  transverse  wake  component  at  d>  =  0.17  which  was  mentioned 
previously.  Since  the  unsteady  pressure  coefficients  Cpand  Cap  are  nondimensionalized 
using  this  transverse  wake  component,  other  effects  such  as  incidence,  airfoil  geometry  and 
impeller  circumferential  pressure  gradients,  which  are  not  considered  in  the  theory,  will 
have  a  greater  relative  effect  on  the  pressure  distribution. 

Figures  6.58  and  6.59  show  the  detuned  cascade  results  for  Nv  =  30,  rj/r0  =  1.10, 
a  =  60°  and  d>  =  0.22.  The  surface  pressure  distributions  show  smoother  behavior  than 
for  the  tuned  cascade,  although  the  theory  does  not  consider  detuning.  The  correlation  of 
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Cap  with  the  theory  is  good  for  both  vanes,  with  the  phase  slightly  higher  for  the  data  than 
the  theory  for  vane  2. 

The  correlation  of  the  data  with  the  theory  is  generally  poor  for  the  tuned  cascade  with 
Nv  =  30.  This  is  believed  to  be  due  to  the  close  spacing  of  the  airfoils,  thus  greatly 
affecting  the  steady  flowfield,  and  distorting  the  wake  which  the  theory  does  not  consider. 
It  is  hypothesized  that  the  detuning  of  the  cascade  lessens  the  effect  of  the  steady  flowfield 
and  correlation  with  the  theory  is  improved. 

Figures  6.60  and  6.61  present  the  unsteady  airfoil  static  pressure  distribution  for 
O  =  0.30  and  0.22  respectively,  with  Nv  =  15,  ri/r0  =  1.10  and  a  =  60°.  The  correlation 
with  theoiy  in  these  cases  is  reasonable  with  some  deviation  in  the  phase  correlation.  The 
theory  generally  appears  to  overpredict  the  pressure  amplitude  from  10  to  25%  chord  for  all 
diffuser  geometries,  with  this  especially  evident  in  the  present  figures.  This  is  possibly 
attributable  to  vane  camber  and  thickness  which  affect  the  steady  flow  around  the  airfoil. 
The  phase  decrease  with  negative  incidence  no  longer  occurs  with  Nv  =  15.  In  Figures 
6.62  and  6.63,  the  same  configuration  with  ri/f0  =  1. 15  is  shown.  The  phase  correlation  is 
good,  with  the  theoretical  amplitude  somewhat  greater  than  the  data.  This  is  attributable  to 
the  unsteady  pressure  being  nonaimensionalized  by  the  wake  transverse  velocity  at 
ri/r0  =  1.10  instead  of  ri/r„  =  1.15,  where  velocity  data  are  not  available. 

Finally,  unsteady  pressure  data  are  presented  for  Nv  =  10  in  Figures  6.64  and  6.65  for 
O  =  0.30  and  0.22  respectively.  The  correlation  with  the  theory  is  poor  in  both  of  these 
figures.  This  behavior  is  not  expected  since  the  correlation  appeared  to  improve  with 
decreasing  vane  number.  The  suction  surface  data  for  Nv  =  10  matches  closely  the  suction 
surface  data  for  Nv  =  15,  however  the  pressure  surface  data  is  considerably  different.  No 
explanation  is  evident  for  this  behavior. 

In  summary,  the  unsteady  lift  increases  with  flow  rate  <t>,  for  Nv  =  30,  with  this 
behavior  particularly  noticed  at  higher  vane  stagger  angles.  Increasing  ri/r0  decreased  the 
unsteady  lift  at  higher  d>  and  a,  although  at  lower  4>,  the  effect  is  less  evident.  With 
Nv  =  15 ,  the  unsteady  lift  is  generally  greater  than  for  Nv  =  30,  with  less  variation  due  to 
flow  rate,  d>.  Increasing  ri/r0  decreased  the  unsteady  lift  significantly.  Nonuniform 
circumferential  spacing  detuning  increased  the  lift  on  one  or  both  of  the  vanes  in  all  cases. 

Poor  correlation  of  the  unsteady  pressure  data  with  the  theory  was  seen  for  Nv  =  30, 
with  the  very  close  vane  spacing  a  probable  cause.  Reasonable  correlation  between  data 
and  theory  was  seen  with  Nv  =  15  at  points  of  low  vane  steady  loading.  An  unexplainable 
poor  correlation  with  theory  was  seen  with  Nv  =  10. 
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6-6  Compressor  Instabilities 

In  all  types  of  pumping  and  compression  systems,  various  forms  of  local  and  system 
instabilities  occur.  Of  primary  relevance  to  gas  compressors  are  surge  and  rotating  stall. 
Surge  is  defined  as  a  global  instability  of  large  mass  flow  oscillation  through  the  machine  at 
low  frequency.  This  can  occur,  as  shown  by  Greitzer  (1981),  when  the  slope  of  the 
compressor  operating  characteristic  (pressure  versus  mass  flow,  as  in  Figures  6.1  through 
6.11)  becomes  positive.  At  this  point,  a  small  perturbation  in  system  mass  flow  can 
increase  in  amplitude  until  the  flow  is  completely  reversed  in  some  systems,  often  causing 
catastrophic  failure. 

Rotating  stall  is  a  more  localized  instability,  in  which  stalling  of  one  or  more  blade 
passages  propagates  around  the  machine  circumference  at  some  fraction  of  the  wheel 
speed,  typically  20  to  70%  according  to  Greitzer  (1981).  Stability  analyses  for  rotating 
stall  are  generally  more  complicated  and  less  exact  than  for  surge,  with  precise  a  priori 
.  prediction  of  a  compressor  stall  point  difficult  to  obtain.  However,  rotating  stall  is 
generally  observed  around  the  zero  slope  portion  of  the  compressor  operating 
characteristic,  with  stalling  more  severe  as  the  flow  rate  is  decreased.  Rotating  stall  is 
sometimes  observed  on  the  negatively  sloped  portion  of  the  compressor  operating 
characteristic,  especially  in  high  speed,  multi-stage  axial  compressors.  Non-stalled 
operation  on  die  positive  portion  of  the  characteristic  is  rare. 

The  classic  explanation  of  rotating  stall  onset  is  given  by  Emmons,  et  al.  (1955).  When 
a  turbomachine  blade  row  is  operating  at  high  incidence  angle,  a  disturbance  in  the  flow  can 
locally  increase  this  angle,  stalling  one  of  the  blades.  As  this  blade  stalls,  the  flow  through 
the  stalled  passage  is  diverted  into  the  adjacent  passages,  decreasing  the  incidence  on  one, 
while  increasing  the  incidence  on  the  other,  thus  causing  it  to  stall.  This  sequence  repeats 
itself,  causing  the  stall  cell  to  rotate  around  the  machine.  More  recent  studies  have  shown 
stall  to  develop  from  small  amplitude  rotating  waves  which  can  increase  in  amplitude  at 
certain  machine  operating  points  .  It  has  been  observed  in  the  literature  (Van  den 
Braembussche,  1984)  and  also  in  the  present  investigation,  that  localized  or  global  blade 
stalling  is  not  always  a  sufficient  condition  for  rotating  stall. 

Surge  and  rotating  stall  are  distinct  phenomena.  However,  they  are  not  unrelated. 
Occurrence  of  rotating  stall  can  often  trigger  the  global  instability  of  surge  in  axial 
compressors.  Although  much  effort  has  been  made  to  predict  and  reduce  these  phenomena 
in  turbomachines,  their  occurrence  and  behavior  is  still  not  well  understood,  especially  in 
centrifugal  compressors.  In  this  regard,  occurrence  of  system  instabilities  in  the  Purdue 
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Research  Centrifugal  Compressor  is  investigated,  with  the  effects  of  vaned  diffuser 
geometry  on  this  behavior  presented. 

A  summary  of  the  various  flow  instabilities  is  given  in  Table  6.2  for  the  vaned  diffuser. 
For  detection  of  surge  and  rotating  stall,  flow  visualization  in  conjunction  with  the  diffuser 
wall  high  response  pressure  transducers  are  utilized.  Flow  visualizing  tufts  are  located  at 
the  impeller  inlet  as  well  as  in  the  radial  diffuser.  Along  with  a  strobe  light,  these  tufts  and 
the  transducer  signals  allow  determination  of  stall  or  surge  occurrence,  location  of  stall 
(impeller  or  diffuser),  relative  amplitude  and  number  of  stall  cells. 


Table  6.2  Compressor  Instabilities 


Ny  =  30 


o=50 0 

0  =  60° 

o  =  70° 

Stall:  ©  =  5.8  Hz  (2  lobes) 
©=12.3  Hz  (5  lobes) 

0  =  0.155 

Surge:  ©  =  5.2  Hz 

0  =  0.148 

Stall:  ©=25  Hz  (1  lobe) 

0  =  0.148 

■ 

Stall:  ©  =  6  Hz  (2  lobes) 
©=  15  Hz  (5  lobes) 

©  =  21  Hz  (3  lobes) 

0  =  0.160 

Surge:  ©  =  5.2  Hz 

0  =  0.148 

No  Stall 

0  =  0.100 

Ny  =  15 


o=50° 

o  =  60° 

& 

11 

& 

Stall:  ©  =  26  Hz  (1  lobe) 

0  =  0.160 

Surge:  ©  =  5.2  Hz 

0  =  0.148 

Stall:  ©  =  25  Hz  (1  lobe) 

O  =  0.145 

Surge:  ©  =  5.2  Hz 

0  =  0.150 

During  compressor  operation  with  the  vaneless  diffuser  (Nv  =  0),  rotating  stall  occurs 
at  0  =  0. 16.  Figure  6.66  shows  the  unsteady  pressure  signal  from  one  of  the  diffuser  wall 
mounted  transducers  as  a  function  of  time.  In  the  bottom  graph,  the  Fast  Fourier 
Transform  (FFT)  of  the  signal  is  shown.  In  the  figure,  one  distinct  frequency  attributable 
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to  rotating  stall  is  seen  at  to  =  26  Hz.  This  frequency  decreases  slightly  as  O  is  reduced. 
The  tufts  at  the  impeller  inlet  and  exit  respond  to  the  stall  cell,  which  appears  to  originate  in 
the  impeller.  At  all  times  during  stall,  only  one  stall  cell  is  observed  (inferred  from  the 
phase  angle  between  the  transducer  signals).  The  machine  is  operated  to  d>  =  0.12,  with 
the  single  cell  impeller  stall  the  only  instability  observed. 

For  vaned  diffuser  operation,  several  types  of  instability  are  observed.  The  frequency 
given  in  Table  2  is  the  cell  passing  frequency,  'ITierefore,  a  stall  pattern  at  to  =  5.8  Hz  with 
2  cells  is  a  2-celled  pattern  rotating  at  2.9  Hz.  The  impeller  rotational  speed  is  29.8  Hz, 
therefore  the  stall  frequency  is  approximately  20%  of  the  impeller  frequency.  When  several 
frequencies  are  present,  some  of  the  harmonics  can  add  or  subtract  to  yield  other 
harmonics. 

The  first  half  of  Table  6.2  show  results  with  Nv  =  30.  In  the  first  row,  flow 
instabilities  with  rj/f0  =  1.10  are  listed.  For  c  =  50°,  the  diffuser  stall  contains  two  distinct 
frequencies  and  number  of  cells.  The  stall  at  to  =  5.8  Hz  appears  to  be  surge  during 
operation.  However,  the  transducer  phase  content  suggests  that  it  is  stall.  Since  rotating 
stall  can  trigger  surge,  it  is  likely  that  this  stall  frequency  will  produce  surge  at  lower  flow 
rates.  For  a  ~  60°,  surge  occurs  at  ®  =  0.148.  The  transducer  signal  is  shown  in  Figure 
6.67.  The  single  surge  frequency  around  to  *  5  Hz  is  seen.  The  phase  difference  between 
the  transducers  is  close  to  zero,  indicating  the  global  mass  flow  oscillation.  No  rotating 
stall  patterns  are  observed  before  surge.  For  a  =  70°,  a  stall  develops  which  is  similar  to 
the  vaneless  diffuser  impeller  stall.  However,  minimal  tuft  movement  is  observed  at  the 
impeller  inlet  as  opposed  to  the  vaneless  case  thus  the  location  of  this  stall  is  uncertain. 

The  behavior  for  rj/r0  =  1.15  is  given  in  the  second  row  of  Table  6.2.  Figure  6.68 
presents  the  transducer  signal  and  FFT  for  o  =  50°.  Three  distinct  frequencies  are  evident 
at  to  =  6,  15  and  21  Hz,  with  the  apparent  number  of  stall  cells  being  2,  5  and  3 
respectively.  It  is  possible  that  2  patterns  exist  and  that  their  sum  or  difference  produces 
the  third  ffequr  .cy.  For  a  =  60°,  the  behavior  is  similar  to  ri/r0  *  1.10,  with  surge 
occurring.  For  a  =  70°,  no  instabilities  are  observed  operating  the  machine  to  d>  =  0.10. 
For  Nv  =  15,  impeller  stall  occurs  at  a  =  50°  and  70°.  The  flow  coefficient  for  a  =  70°  is 
lower,  d>  =  0.145,  than  for  the  vaneless  impeller  stall  where  =  0.160.  For  a  =  60°, 
surge  occurs  as  in  the  previous  cases. 

The  vaned  diffuser  is  seen  to  have  a  significant  influence  on  the  occurrence  of 
compressor  instabilities.  The  flow  coefficient  for  stall  can  be  reduced  by  proper  geometry 
selection.  For  Nv  =  30,  o  =  70°  and  ri/rQ  =  1.15,  stall  and  surge  behavior  are  completely 
eliminated,  thus  giving  support  to  the  belief  that  rotating  disturbances  will  only  be  amplified 
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under  certain  conditions,  with  static  stalling  of  individual  blades  or  vanes  not  always 
sufficient 


Pressure  Rise  Coefficient, 


Flow  Coefficient,  O 

Centrifugal  Compressor  Performance  with  Vaneless  Diffuser 


Figure  6.1 
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Vane  Number  on  Compressor  Performance  (ri/r0  = 


Figure  6.7  Effect  of 
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Figure  6. 1 1  Detuned 


Impeller  Radius  Radio,  r/r0 


Figure  6.12  Impeller  Static  Pressure  (Nv  =  0,  d>  =  0.17) 


Figure  6.14 


115 


Cap 


n - 1 - 1 

0.6  0.8  1.0 


Impeller  Radius  Radio,  r/r0 


r  Static  Pressure  (Nv  =  0,0  =  0.30) 


Impeller  Blade  Surface  Static  Pressure, 
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Impeller  Radius  Radio,  r/r0 

Figure  6.15  Impeller  Static  Pressure  (rj/r0  =  1.10,  a  =  50°,  Nv  =  30, 0  =  0.30) 


Impeller  Blade  Surface  Static  Pressure,  Cp  and  C^p 


Impeller  Radius  Radio,  r/r0 


Figure  6.16  Impeller  Static  Pressure  (ri/r0  =  1.10,  c  =  60°,  Nv  =  30,  d>  =  0.22) 


118 


Figure  6.17  Incidence  Effects  on  Diffuser  Vane  Steady  Static  Pressure 
(ri/r0  =  1.10,  a  =  60°,  Nv  =  15) 
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Figure  6.18  Vane  Number  Effects  on  Diffuser  Vane  Steady  Static  Pressure 
(ri/r0  =  1.10,0  =  60°,  <D  =  0.22) 
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Figure  6.19  Detuned  Diffuser  Vane  Steady  Static  Pressure 
(ri/r0  =  1.10,  c  =  60°,  Nv  =  30, 0  =  0.30) 


Figure  6.21  LDV  Circumferential  Measurement  Bins 


Steady  Tangential  Velocity,  iyU6I  Steady  Radial  Velocity,  U/Url 
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Figure  6.23  Vaneless  Diffuser  Time-Averaged  Velocity  (d>  =  0.30) 


125 


Figure  6.24 
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Vaneless  Diffuser  Time- Averaged  Velocity  (O  =  0.22) 
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Axial  Position  =  Tip 


Figure  6.26  Vaneless  Diffuser  Total  Velocity  (Tip  Position,  O  =  0.30) 
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Figure  6.29  Vaneless  Diffuser  Total  Velocity  (Tip  Position,  <D  =  0.22) 
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Figure  6.33  Vaneless  Diffuser  Total  Velocity  (Mean  Position,  0  =  0.17) 
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Figure  6.34  Vanelcss  Diffuser  Total  Velocity  (Hub  Position,  O  =  0.17) 
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Figure  6.35  Vaneless  Diffuser  Unsteady  Velocity  (Tip  Position,  O  =  0.30) 
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Figure  6.36  Vaneless  Diffuser  Unsteady  Velocity  (Mean  Position,  O  =  0.30) 
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Figure  6.37  Vaneless  Diffuser  Unsteady  Velocity  (Hub  Position,  O  =  0.30) 
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Figure  6.38  Vaneless  Diffuser  Unsteady  Velocity  (Tip  Position,  <D  =  0.22) 
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Figure  6.42  Vaneless  Diffuser  Unsteady  Velocity  (Mean  Position,  0  =  0. 17) 
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Figure  6.43  Vaneless  Diffuser  Unsteady  Velocity  (Hub  Position,  <D  =  0.17) 
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Figure  6.46  Vaneless  Diffuser  First  Harmonic  Unsteady  Velocity  and  Theory  (<t>  =  0.22) 
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Figure  6.48  Vaneless  Diffuser  First  Harmonic  Unsteady  Velocity  and  Theory 
(Theory  <j)w  =  -160°,O  =  0.30) 
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Figure  6.52  Diffuser  Vane  First  Harmonic  Unsteady  Lift  (rj/r0  =  1. 15,  Nv 
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Figure  6.55  Diffuser  Vane  First  Harmonic  Unsteady  Static  Pressure  (ri/rD  =  1.10,  Nv  =  30,  a  =  60°,  <E>  =  0.30) 


Figure  6.56  Diffuser  Vane  First  Harmonic  Unsteady  Static  Pressure  (ri/rQ  =  1.10,  Nv  =  30,  a  =  60°,  0  =  0.22) 


Figure  6.57  Diffuser  Vane  First  Harmonic  Unsteady  Static  Pressure  (ri/r0  =  1.10,  Nv  =  30,  a  =  60°,  0  =  0.17) 
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Figure  6.58  Diffuser  Vane  First  Harmonic  Unsteady  Static  Pressure  (Vane  1, 40%  detM  x\fr0  =  1.10,  Nv  =  30,  o  =  60°,  O  =  0.22) 


Diffuser  Vane  First  Harmonic  Unsteady  Static  Pressure  (Vane  2, 40%  del.,  n/r0  =  1.10,  Nv  =  30,  c  =  60°,  0  =  0.22) 


Diffuser  Vane  First  Harmonic  Unsteady  Static  Pressure  (ri/rc  =  1.10,  Nv  =  15,  o  =  60°,  O  =  0.30) 
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Figure  6.62  Diffuser  Vane  First  Harmonic  Unsteady  Static  Pressure  (ri/r0  =  1.15,  Nv  =  15,  a  =  60°,  O  =  0.30) 
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Diffuser  Vane  First  Harmonic  Unsteady  Static  Pressure  (rj/r0  =  1.15,  Nv  =  15,  c  =  60°,  <I>  =  0.22) 
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Figure  6.64  Diffuser  Vane  First  Harmonic  Unsteady  Static  Pressure  (ri/r0  =  1.10,  Nv  =  10,  a  =  60°,  O  =  0.30) 
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Figure  6.67  Rotating  Stall  Unsteady  Static  Pressure 

(rj/r0  =  1.10,  Nv  =  30, a=  60°  <D  =  0.148) 


CHAPTER  7 

SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


An  investigation  of  some  of  the  unsteady  phenomena  occurring  in  centrifugal 
compressors  has  been  performed,  with  the  diffuser  vane  unsteady  pressure  generated  by 
impeller  wakes  of  particular  interest  It  was  found  that  little  attention  has  been  devoted  in 
the  literature  to  the  effect  of  the  vaned  diffuser  on  the  unsteady  aerodynamics  inherent  in 
centrifugal  compressors.  Therefore,  analyses  were  developed  to  predict  the  wake  behavior 
and  diffuser  vane  unsteady  loading.  In  conjunction  with  these,  a  series  of  experiments  was 
performed  to  evaluate  this  unsteady  behavior.  The  effect  of  vaned  diffuser  geometry  on  the 
compressor  steady  performance  and  loading  as  well  as  surge  and  rotating  stall  was 
experimentally  determined. 

The  performance  of  the  Purdue  Research  Centrifugal  Compressor  was  measured  for 
various  diffuser  geometries.  It  was  found  that  the  vaned  duruser  increases  efficiency  and 
overall  pressure  rise  with  proper  geometry  selection.  The  diffuser  stagger  angle  and 
number  of  diffuser  vanes  were  seen  to  be  particularly  important,  with  a  smaller  number  of 
diffuser  vanes  generally  preferable.  The  mechanism  of  this  performance  improvement  was 
investigated  by  considering  the  rotating  impeller  blade  and  diffuser  vane  state  loading.  The 
vaned  diffuser  was  determined  to  affect  the  impeller  loading  only  at  higher  flow  rates  with 
small  vane  loading.  The  diffuser  loading  increased  with  steady  incidence  and  with  a 
decreasing  number  of  vanes.  It  was  concluded  that  the  vaned  diffuser  affects  the  aft 
impeller  and  vaneless  space  flowfields,  reducing  mixing  and  inhibiting  secondary  flows, 
thereby  improving  performance.  Performance  improvements  were  only  realized  when  the 
diffuser  vanes  were  lightly  loaded,  with  no  performance  advantage  from  a  loaded  diffuser. 

The  wake  velocities  in  the  vaneless  radial  diffuser  were  experimentally  determined 
using  laser  anemometry.  The  radial  and  tangential  unsteady  wake  velocities  decreased 
rapidly  with  increasing  radius,  with  a  greater  decrease  at  lower  flow  rates.  A  theoretical 
model  to  predict  this  behavior  was  developed  by  obtaining  a  solution  to  the  linearized  Euler 
Equations  in  a  radial  vaneless  diffuser  for  converted  rotational  wakes.  At  higher  flow 
rates,  the  theory  matched  the  data  for  the  radial  velocity  in  a  general  sense. 
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However,  the  measured  tangential  velocity  decreased  more  rapidly  with  increasing 
radius  than  predicted  by  the  theory.  The  amplitude  ratio  and  phase  angle  between  the  wake 
radial  and  tangential  velocities  was  concluded  to  affect  the  correlation  with  theory,  as  well 
as  an  unsteady  circumferential  static  pressure  generated  by  the  rotating  impeller. 

Unsteady  surface  static  pressure  measurements  were  made  on  the  diffuser  vanes  to 
determine  unsteady  loading.  The  unsteady  lift  increased  with  flow  rate  using  30  diffuser 
vanes,  especially  at  higher  vane  stagger.  Increasing  the  diffuser  leading  edge/impeller  exit 
radial  spacing  decreased  the  lift  at  higher  flow  rates  and  stagger  angles,  although  at  lower 
flow  rates,  the  effect  was  less  pronounced.  With  15  diffuser  vanes,  the  unsteady  lift  was 
generally  greater,  with  less  variation  due  to  flow  rate.  Increasing  the  diffuser  leading 
edge/impeller  exit  radial  spacing  dec  rented  the  unsteady  lift  significantly.  Nonuniform 
circumferential  spacing  detuning  increased  the  lift  on  one  or  both  of  the  vanes  in  all  cases. 

A  theory  to  predict  the  unsteady  diffuser  vane  loading  was  developed  by  applying 
conformal  mapping  »o  existing  axial  cascade  unsteady  aerodynamic  theory  and 
incorporating  the  wake  model  as  a  boundary  condition.  Poor  correlation  of  the  data  with 
the  theory  was  seen  with  30  diffuser  vanes,  which  was  attributed  to  the  very  close  vane 
spacing.  Reasonable  correlation  between  data  and  theory  was  seen  with  15  diffuser  vanes 
at  points  of  low  vane  steady  loading.  An  unexplainable  poor  correlation  with  theory  was 
seen  using  10  diffuser  vanes.  Mismatch  of  data  with  theory  is  believed  to  be  due  to 
mismatch  of  the  wake  model  and  steady  loading  on  the  diffuser  vanes,  which  the  theory 
does  not  consider. 

Finally,  flow  instabilities  in  the  compressor  were  experimentally  investigated.  A  one- 
Iobed  impeller  stall  existed  at  low  flow  coefficient  for  the  vaneless  diffuser.  For  the  vaned 
diffuser,  several  instabilities  were  observed.  With  30  diffuser  vanes,  diffuser  stall  and 
surge  were  the  only  instabilities  occurring.  For  60°  vane  stagger,  surge  typically  occurred 
at  a  flow  coefficient  of  approximately  0. 155,  with  rotating  stall  not  observed.  At  the  higher 
diffuser  vane  leading  edge/impeller  exit  ratio  for  60°  vane  stagger,  the  compressor  was 
operated  at  very  low  flow  coefficient  without  surge  or  rotating  stall.  With  15  diffuser 
vanes,  the  impeller  stalled  before  the  diffuser,  with  surge  again  occurring  at  60°  vane 
stagger  without  rotating  stall  occurrence..  It  was  seen  that  localized  vane  stalling  was  not  a 
sufficient  condition  for  rotating  stall  to  occur.  Also,  significant  improvement  in  stall 
margin  was  realized  by  proper  adjustment  of  the  vaned  diffuser. 

The  recommendations  for  further  study  are  basically  twofold.  First,  further 
experimental  investigation  is  needed  to  more  accurately  determine  the  parameters  that  affect 
the  impel ler/radial  diffuser  unsteady  wake  behavior.  The  determination  of  how  the 
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amplitude  ratio  and  phase  angle  between  the  wake  velocities  affect  the  wake  behavior  and 
how  these  vary  with  impeller  geometry,  operating  point  and  roi  ional  speed  would  be 
useful.  More  data  on  unsteady  diffuser  vane  loading  would  also  be  beneficial  to  more 
accurately  determine  the  effects  of  diffuser  vane  row  geometry  on  the  unsteady  response. 
More  evidence  on  the  effect  of  diffuser  geometry  on  rotating  stall  would  also  be  a 
significant  experimental  contribution.  Since  the  current  investigation  was  performed  at  low 
speed,  it  is  important  to  establish  the  correlation  between  measurements  made  at  low  speed 
and  behavior  of  high  speed  machines. 

The  second  recommendation  concerns  the  theoretical  prediction  of  unsteady  behavior. 
A  more  accurate  model  for  radial  diffuser  wake  behavior  would  certainly  be  useful.  It  is 
likely  a  model  that  considers  the  circumferential  pressure  variation  as  a  boundary  condition 
will  be  necessary.  Once  this  is  accomplished,  incorporating  this  model  into  unsteady  vane 
loading  prediction  is  of  primary  importance.  The  effects  of  vane  loading  and  gust 
distortion  should  be  considered  in  the  unsteady  loading  analysis.  Since  the  unsteady 
diffuser  vane  loading  analysis  is  only  valid  for  incompressible  flow,  prediction  of  unsteady 
loading  in  high  speed  machines  will  require  consideration  of  compressibility  effects. 

In  conclusion,  a  more  significant  data  base  along  with  more  general  analyses  will  be 
necessary  before  accurate  prediction  of  unsteady  aerodynamics  in  centrifugal  compressors 
is  realized. 
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